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Set Name Query Hit Count Set Name 

side by side result set 

DB=USPT,PGPBJPAB,EPAB,DWPI,TDBD; PLUR=YES; OP=ADJ 

L 19 L 1 7 and ((optim$9 or ideal$5 or best) with (imag$6)) 84 LI 9 

L18 L17 and ((optim$9 or ideal$5 or best) with (imag46)) 0 L18 

LI 6 and (boundar$6 or artifact$3 or artefact$3 or ghost$5 or blur$6 ^ ^7 
or alias$6 or slab or edg$4) 

L14 and ("set" or group$5 or plurality or sub-set or subset or "sub ^ ^jg 
set") 

L15 L14 and ("set" or group or plurality or sub-set or subset or "sub set") 104 L15 

LI 3 and (kspace or k-space or "k space" or "kx" or "ky" or "kz" or ^ ^4 
raw or ((image or frequency) with space) ) 

L13 L12and(optim$9orideal$5orbest) 121 U3 

T , 0 LI 1 and ((complet$5 or entir46 or total$3 or finish$5 or whole) with 

^ (imag$6)) 144 ^ 

^ j L10 and ((complet$5 or entir46 or total$3 or finish$5 or whole) with ^ 4g ^ 
((field of view) or field-of-view or fov)) 

L10 L9 and (complet$5 or entir46 or total$3 or finish$5 or whole) 526 LI 0 

L9 L8 and ((imag$6) with ((field of view) or field-of-view or fov)) 556 L9 

L8 L7 and (comput$6 or proces$9 or procces$9 or program$8) 840 L8 

L7 L6 and ((imag$6) with (scan$7 or apparatus or device)) 840 L7 

L6 L5 and (comput$6 or proces$9 or procces49 or program$8) 944 L6 

L5 L4 and (scan$7 or apparatus or device) 974 L5 

L3 and (direction$3 or axes or axis or encod$6 or increment$5 or 

L4 step$g) 1U1U LA 

L3 L2 and (imag$6) 1028 L3 

L2 LI and (gradient) 1069 L2 

LI ((field of view) or field-of-view or fov) 5084 LI 
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Search Results - Record(s) 1 through 84 of 84 returned. 
□ 1. Document ID: US 20030189098 Al 

L19: Entry 1 of 84 File: PGPB Oct 9, 2003 



PGPUB- DOCUMENT -NUMBER : 20030189098 
PGPUB- FILING -TYPE : new 

DOCUMENT- IDENTIFIER: US 20030189098 Al 

TITLE: Planar light illumination and imaging (PLUM) systems employing LED-based 
planar light illumination arrays (PLIAS) and linear electronic image detection 
arrays 

PUBLICATION-DATE: October 9, 2003 

INVENTOR- INFORMATION : 
NAME 

Tsikos, Constantine J. 
Good, Timothy A. 
Schnee, Michael D. 
Zhu, Xiaoxun 
Amundsen, Thomas 
Knowles, C. Harry 

US - CL - CURRENT : 235/454 



CITY 

Voorhees 

Clementon 

Aston 

Marl ton 

Turnersville 

Moorestown 



STATE COUNTRY RULE-47 

NJ US 

NJ US 

PA US 

NJ US 

NJ US 

NJ US 




! 2. Document ID: US 20030150917 Al 

L19: Entry 2 of 84 File: PGPB Aug 14, 2003 



PGPUB -DOCUMENT -NUMBER: 20030150917 
PGPUB -FILING -TYPE : new 

DOCUMENT- IDENTIFIER: US 20030150917 Al 



TITLE: Planar light illumination and imaging (PLUM) system employing led-based 
planar light illumination arrays (PLIAS) and an area-type image detection array 

PUBLICATION-DATE: August 14, 2003 



INVENTOR- INFORMATION : 
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RecoAl List Display 



http://westhrs:800^in/gate.exe?5=TOG 



PT,PGPB,JPAB > EPAB > DWPI,TDBD&ESN T AME^ 



NAME 


CITY 


STATE 


COUNTRY 


Tsikos, Constantine J. 


Voorhees 


NJ 


US 


Good/ Timothy 


Clementon 


NJ 


US 


Schnee, Michael D. 


Aston 


PA 


US 


Zhu, Xiaoxun 


Marlton 


NJ 


US 


Amundsen, Thomas 


Turnersville 


NJ 


US 


Knowles, C. Harry 


Moorestown 


NJ 


US 



RULE-47 



US - CL - CURRENT : 23 5/454 



1 



O 3. Document ID: US 20030150916 Al 

L19: Entry 3 of 84 File: PGPB 



Aug 14, 2003 



PG PUB - DOCUMENT - NUMBER : 20030150916 
PGPUB- FILING -TYPE : new 

DOCUMENT- IDENTIFIER: US 20030150916 Al 



TITLE: LED-based planar 


light 


illumination and 


imaging 


(PLUM) systems 


PUBLICATION-DATE: August 


14, 


2003 






INVENTOR- INFORMATION : 










NAME 




CITY 


STATE COUNTRY RULE - 4 7 


Tsikos, Constantine J. 




Voorhees 


NJ 


US 


Good, Timothy 




Clementon 


NJ 


US 


Schnee, Michael D. 




Aston 


PA 


US 


Zhu, Xiaoxun 




Marlton 


NJ 


US 


Amundsen, Thomas 




Turnersville 


NJ 


US 


Knowles, C. Harry 




Moorestown 


NJ 


us 


US - CL - CURRENT : 235/454 
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O 4. Document ID: US 20030146282 Al 

L19: Entry 4 of 84 File: PGPB Aug 7, 2003 

PG PUB - DO CUMENT - NUMBE R : 20030146282 
PGPUB- FILING -TYPE : new 

DOCUMENT- IDENTIFIER: US 2003 0146282 Al 

TITLE: Bioptical product and produce identification systems employing planar laser 
illumination and imaging (PLIM) based subsystems 

PUBLICATION-DATE: August 7, 2003 

I NVENTOR - I NFORMAT I ON : 
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NAME 


CITY 


STATE 


COUNTRY 


Tsikos, Constant ine J. 


Voorhees 


NJ 


US 


Wirth, Allan 


Bedford 


MA 


US 


Good, Timothy A. 


Clementon 


NJ 


US 


Jankevics, Andrew 


Westford 


MA 


US 


Kim, Steve Y. 


Cambridge 


MA 


US 


Amundsen, Thomas 


Turnersville 


NJ 


us 


Naylor, Charles A. 


Sewell 


NJ 


us 


Dobbs, Russell Joseph 


Cherry Hill 


NJ 


us 


Giordano, Patrick A. 


Blackwood 


NJ 


us 


Yorsz, Jeffery 


Winchester 


MA 


us 


Schmidt, Mark S. 


Williamstown 


NJ 


us 


Colavito, Stephen J. 


Brookhaven 


PA 


us 


Wilz, David M. SR. 


Sewell 


NJ 


us 


Au, Ka Man 


Philadelphia 


PA 


us 


Svedas, William 


Deptford 


NJ 


us 


Ghosh, Sankar 


Glenolden 


PA 


us 


Schnee, Michael D . 


Aston 


PA 


us 


Zhu, Xiaoxun 


Marl ton 


NJ 


us 


Knowles , C . Harry 


Moor es town 


NJ 


us 



RULE-47 



US - CL - CURRENT : 235/454 
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O 5. DocumentID: US 20030135111 Al 

L19: Entry 5 of 84 File: PGPB 



Jul 17, 2003 



PG PUB -DOCUMENT -NUMBER: 20030135111 
PG PUB-FILING-TYPE : new 

DOCUMENT- IDENTIFIER: US 20030135111 Al 

TITLE: Method and apparatus for magnetic resonance arteriography using contrast 
agents 

PUBLICATION-DATE: July 17, 2003 

INVENTOR- INFORMATION : 
NAME 

Meaney, James F.M. 
Prince, Martin R. 



CITY 
Leeds 
Ann Arbor 



STATE COUNTRY RULE-47 

MI GB 
US 



US - CL - CURRENT : 600/422 
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O 6. DocumentID: US 20030114748 Al 

L19: Entry 6 of 84 File: PGPB 



Jun 19, 2003 
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PG PUB -DOCUMENT -NUMBER: 2003 0114748 
PGPUB- FILING -TYPE : new 

DOCUMENT- IDENTIFIER: US 20030114748 Al 

TITLE: Inherently de-coupled sandwiched solenoidal array coil 
PUBLICATION-DATE: June 19, 2003 



INVENTOR- INFORMATION : 
NAME 

Su, Sunyu 
Kaufman, Leon 
Arakawa, Mitsuaki 
Carlson, Joseph W. 

US - CL - CURRENT : 600/422 



CITY 

South San Francisco 
San Francisco 
Hillsborough 
Kensington 



STATE 

CA 

CA 

CA 

CA 



COUNTRY 

US 

US 

US 

US 



RULE-47 
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O 7. Document ID: US 20030109782 Al 

L19: Entry 7 of 84 File: PGPB 

PG PUB - DOCUMENT - NUMBER : 2003 0109782 
PGPUB- FILING- TYPE : new 

DOCUMENT -IDENTIFIER: US 20030109782 Al 

TITLE: Inherently de-coupled sandwiched solenoidal array coil 
PUBLICATION-DATE: June 12, 2003 
INVENTOR- INFORMATION : 



NAME 

Su, Sunyu 
Kaufman, Leon 



CITY 

South San Francisco 
San Francisco 



STATE COUNTRY 
CA US 
CA US 



Jun 12, 2003 



RULE-47 



US - CL - CURRENT : 6 00/421 
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8. Document ID: US 20030102379 Al 

L19: Entry 8 of 84 File: PGPB 



Jun 5, 2003 



PG PUB - DOCUMENT - NUMBER : 20030102379 
PGPUB- FILING- TYPE : new 

DOCUMENT- IDENTIFIER: US 20030102379 Al 

TITLE: LED-based planar light illumination and 

PUBLICATION-DATE: June 5, 2 003 



imaging (PLUM) engine 
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I NVENTOR - I NFORMAT I ON : 



NAME* 


r 1 T TV 


CTATD 


L.UUJN IKi 


Tsikos , Constant ine J . 


Voorhees 


VT T 
JNI J 


T TO 

Ub 


Knowles / C . Harry 


Moor est own 


"NT T 

vi a 


T TO 

Ub 


Good, Timothy A. 


Clementon 




T TO 

Ub 


Schnee, Michael D. 


Aston 


PA 


US 


Zhu, Xiaoxun 


Marlton 


NJ 


us 


Amundsen, Thomas 


Turnersville 


NJ 


us 


Schmidt, Mark C. 


Williamstown 


NJ 


us 


Giordano, Patrick A. 


Blackwood 


NJ 


us 



RULE -47 



US - CL - CURRENT : 235/462 ^45 
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O 9. Document ID: US 20030100825 Al 

L19: Entry 9 of 84 File: PGPB 

PG PUB - DO CUMENT - NUMBER : 2003 0100825 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20030100825 Al 



TITLE: METHOD AND SYSTEM FOR EXTENDED VOLUME 
PUBLICATION-DATE: May 29, 2003 



IMAGING USING MR I 



INVENTOR- INFORMATION : 
NAME 

Demoulin, Charles Lucian 
Zhu , Yudong 

US -CL- CURRENT: 6 00/410 



CITY STATE COUNTRY 

Ballston Lake NY US 

Clifton Park NY US 



May 29, 2003 



RULE- 4 7 



m£ge : j 



May 29, 2003 



□ 10. Document ID: US 20030098353 Al 

L19: Entry 10 of 84 File: PGPB 

PGPUB- DOCUMENT -NUMBER: 20030098353 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20030098353 Al 

TITLE: Planar laser illumination and imaging (PLUM) engine 
PUBLICATION-DATE: May 29, 2003 
INVENTOR- INFORMATION : 
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NAME 


CITY 


STATE 


COUNTRY 


Tsikos, Constantine J. 


Voorhees 


NJ 


US 


Knowles , C . Harry 


Moor es town 


NJ 


US 


Zhu, Xiaoxun 


Marlton 


NJ 


US 


Schnee, Michael D. 


Aston 


PA 


US 


Amu nd s e n , Thoma s 


Turnersville 


NJ 


US 


Schmidt, Mark C. 


Williamstown 


NJ 


US 


Giordano, Patrick A. 


Blackwood 


NJ 


US 



RULE -4 7 



US - CL - CURRENT : 235/472.01 
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O 11. Document ID: US 20030098349 Al 

L19: Entry 11 of 84 File: PGPB 



May 29, 2003 



PG PUB -DOCUMENT -NUMBER: 20030098349 
PGPUB- FILING -TYPE : new 

DOCUMENT- IDENTIFIER: US 20030098349 Al 



imaging (PLUM) devices 



TITLE: Hand -supportable planar laser illumination and 
employing linear electronic image detection arrays and optically-combined planar 
laser illumination beams (PLIBS) produced from a multiplicity of laser diode sources 
to achieve a reduction in speckle-pattern noise power in said PLUM devices 

PUBLICATION-DATE: May 29, 2003 



I NVENTOR - I NFORMAT I ON : 
NAME 

Tsikos, Constantine J. 
Knowles, C. Harry 
Zhu, Xiaoxun 
Schnee, Michael D. 
Good, Timothy A. 



CITY 


STATE 


COUNTRY 


Voorhees 


NJ 


US 


Moorestown 


NJ 


US 


Marlton 


NJ 


US 


Aston 


PA 


US 


Clement on 


NJ 


US 



RULE -4 7 



US - CL - CURRENT : 235 /462.01 
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O 12. Document ID: US 20030094495 Al 

L19: Entry 12 of 84 File: PGPB May 22, 2003 

PGPUB -DOCUMENT -NUMBER : 2003 0094495 
PGPUB-FILING-TYPE : new 

DOCUMENT- IDENTIFIER: US 20030094495 Al 

TITLE: Nuclear resonance based scanning system having an automatic object 
identification and attribute information acquisition and linking mechanism 
integrated therein 
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PUBLICATION-DATE: May 22, 2003 



I NVENTOR - I NFORMAT I ON : 
NAME 

Knowles, C. Harry 
Schmidt, Mark C. 
Fisher, Dale 



CITY 

Moor es town 

Williamstown 

Voorhees 



STATE 
NJ 
NJ 
NJ 



COUNTRY 

US 

US 

US 



RULE -47 



US -CL- CURRENT: 235/462. 14 
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sequences; 



...J 13. Document ID: US 20030089778 Al 

L19: Entry 13 of 84 File: PGPB 



May 15, 2003 



Images of a moving 



PG PUB -DOCUMENT -NUMBER: 2003 0089778 
PGPUB- FILING -TYPE : new 

DOCUMENT- IDENTIFIER: US 20030089778 Al 

TITLE: Method of and system for automatically producing digital 
object, with pixels having a substantially uniform white level independent of the 
velocity of said moving object 

PUBLICATION-DATE: May 15, 2003 

INVENTOR- INFORMATION : 

RULE- 4 7 



NAME 


CITY 


STATE 


COUNTRY 


Tsikos, Constantine J. 


Voorhees 


NJ 


US 


Knowles, C. Harry 


Moorestown 


NJ 


US 


Zhu, Xiaoxun 


Marlton 


NJ 


US 


Schnee, Michael D. 


Aston 


PA 


US 


Au, Ka Man 


Philadelphia 


PA 


US 


Wirth, Allan 


Bedford 


MA 


US 


Good, Timothy A. 


Clementon 


NJ 


US 


Ghosh, Sankar 


Glenolden 


PA 


US 


Amundsen, Thomas 


Turnersville 


NJ 


US 



US - CL - CURRENT : 235/454 
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□ 14. Document ID: US 20030085281 Al 

L19: Entry 14 of 84 File: PGPB May 8, 2003 

PG PUB-DOCUMENT-NUMBER : 2003 0085281 
PGPUB -FILING -TYPE : new 

DOCUMENT- IDENTIFIER: US 20030085281 Al 

TITLE: Tunnel-type package identification system having a remote image keying 
station with an ethernet -over- fiber-optic data communication link 
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PUBLICATION-DATE: May 8, 2003 

INVENTOR- INFORMATION : 

NAME CITY 

Knowles, C. Harry Moorestown 

Kim, Steven Y. Cambridge 



STATE 

NJ 

MA 



COUNTRY 

US 

US 



RULE- 47 



US -CL- CURRENT: 235/454 
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O 15. Document ID: US 20030085280 Al 

L19: Entry 15 of 84 File: PGPB 



May 8, 2003 



PG PUB -DOCUMENT -NUMBER: 20030085280 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20030085280 Al 

TITLE: Method of speckle-noise pattern reduction and apparatus therefor based on 
reducing the spatial -coherence of the planar laser illumination beam before it 
illuminates the target object by applying spatial intensity modulation techniques 
during the transmission of the PLIB towards the target 

PUBLICATION-DATE: May 8, 2 003 

I NVENTOR - I NFORMAT I ON : 

RULE-47 



NAME 


CITY 


STATE 


COUNTRY 


Tsikos, Constantine J. 


Voorhees 


NJ 


US 


Knowles , C . Harry 


Moorestown 


NJ 


US 


Wirth, Allan 


Bedford 


MA 


US 


Good, Timothy A. 


Clementon 


NJ 


US 


Jankevics, Andrew 


Westford 


MA 


US 



US - CL - CURRENT : 235 / 454 



O 16. Document ID: US 20030080192 Al 

L19: Entry 16 of 84 File: PGPB May 1, 2003 

PG PUB -DOCUMENT -NUMBER: 20030080192 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20030080192 Al 

TITLE: Neutron-beam based scanning system having an automatic object identification 
and attribute information acquisition and linking mechanism integrated therein 

PUBLICATION-DATE: May 1, 2003 

I NVENTOR - I NFORMAT I ON : 
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lUSPT,PGPB,JPABJEPABJ)WPI,TDBD&ESNAME=- 



NAME 


CITY 


STATE 


COUNTRY 


Tsikos, Constantine J. 


Voorhees 


NJ 


US 


Knowles , C . Harry 


Moorestown 


NJ 


US 


Zhu, Xiaoxun 


Marlton 


NJ 


US 


Schnee, Michael D. 


Aston 


PA 


US 


Au, Ka Man 


Philadelphia 


PA 


US 


Wirth, Allan 


Bedford 


MA 


US 


Good, Timothy A. 


Clementon 


NJ 


US 


Jankevics, Andrew 


Westf ord 


MA 


US 


Ghosh, Sankar 


Glenolden 


PA 


US 


Naylor, Charles A. 


Sewell 


NJ 


US 


Amundsen, Thomas 


Turnersville 


NJ 


US 


Blake, Robert 


Woodbury Heights 


NJ 


US 


Svedas, William 


Deptford 


NJ 


US 


Defoney, Shawn 


Runnemede 


NJ 


US 


Skypala, Edward 


Blackwood 


NJ 


us 


Vatan, Pirooz 


Wilmington 


MA 


us 


Dobbs, Russell Joseph 


Cherry Hill 


NJ 


us 


Kolis, George 


Pennsauken 


NJ 


us 


Schmidt, Mark C. 


Williamstown 


NJ 


us 


Yorsz, Jeffery 


Winchester 


MA 


us 


Giordano, Patrick A. 


Blackwood 


NJ 


us 


Colavito, Stephen J. 


Brookhaven 


PA 


us 


Wilz, David W. SR. 


Sewell 


NJ 


us 


Schwartz, Barry E. 


Haddonf ield 


NJ 


us 


Kim, Steven Y. 


Cambridge 


MA 


us 


Fisher, Dale 


Voorhees 


NJ 


us 


Tassell, Jon Van 


Winchester 


MA 


us 



RULE-47 



US - CL - CURRENT : 235/462 . 14 
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□ 17. Document ID: US 20030080190 Al 

L19: Entry 17 of 84 File: PGPB May 1, 2003 

PGPUB - DOCUMENT - NUMBER : 20030080190 
PGPUB- FILING-TYPE : new 

DOCUMENT- IDENTIFIER: US 20030080190 Al 

TITLE: Method of and system for automatically producing digital images of a moving 
object, with pixels having a substantially uniform white level independent of the 
velocity of said moving object 

PUBLICATION-DATE: May 1, 2003 

I NVENTOR - I NFORMAT I ON : 
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NAME 


CITY 


STATE 


COUNTRY 


Tsikos, Constantine J. 


Voorhees 


NJ 


US 


Knowles , C . Harry 


Moore s town 


NJ 


US 


Zhu, Xiaoxun 


Marl ton 


NJ 


US 


Schnee, Michael D. 


Aston 


PA 


us 


Au, Ka Man 


Philadelphia 


PA 


us 


Wirth, Allan 


Bedford 


MA 


us 


Good, Timothy A. 


Clementon 


NJ 


us 


Jankevics, Andrew 


Westf ord 


MA 


us 


Ghosh, Sankar 


Glenolden 


PA 


us 


Naylor, Charles A. 


Sewell 


NJ 


us 


Amundsen, Thomas 


Turnersville 


NJ 


us 


Blake, Robert 


Woodbury Heights 


NJ 


us 


Svedas, William 


Deptford 


NJ 


us 


Defoney, Shawn 


Runnemede 


NJ 


us 


Skypala, Edward 


Blackwood 


NJ 


us 


Vatan, Pirooz 


Wilmington 


MA 


us 


Dobbs, Russell Joseph 


Cherry Hill 


NJ 


us 


Kolis, George 


- Pennsauken 


NJ 


us 


Schmidt, Mark C. 


Williamstown 


NJ 


us 


Yorsz, Jeffery 


Winchester 


MA 


us 


Giordano, Patrick A. 


Blackwood 


NJ 


us 


Colavito, Stephen J. 


Brookhaven 


PA 


us 


Wilz, David W. SR. 


Sewell 


NJ 


us 


Schwartz, Barry E. 


Haddonf ield 


NJ 


us 


Kim, Steven Y. 


Cambridge 


MA 


us 


Fisher, Dale 


Voorhees 


NJ 


us 


Tassell, Jon Van 


Winchester 


MA 


us 



RULE-47 



US -CL- CURRENT: 235/462 .03 




O 18. Document ID: US 20030071 128 Al 

L19: Entry 18 of 84 File: PGPB Apr 17, 2003 

PG PUB - DOCUMENT - NUMBER : 20030071128 
PGPUB- FILING -TYPE ; new 

DOCUMENT- IDENTIFIER: US 2003 0071128 Al 

TITLE: Led-based planar light illumination and imaging (PLUM) based camera system 
employing real-time object coordinate acquistion and producing to control automatic 
zoom and focus imaging optics 

PUBLICATION-DATE: April 17, 2003 

INVENTOR- INFORMATION : 
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NAME 


CITY 


STATE 


COUNTRY 


Tsikos, Constantine J. 


Voorhees 


NJ 


US 


Knowles , C . Harry 


Moorestown 


NJ 


US 


Zhu, Xiaoxun 


Marl ton 


NJ 


US 


Schnee, Michael D. 


Aston 


PA 


US 


Au, Ka Man 


Philadelphia 


PA 


US 


Good, Timothy A. 


Clement on 


NJ 


US 


Ghosh, Sankar 


Glenolden 


PA 


US 


Amundsen, Thomas 


Turnersville 


NJ 


US 
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US -CL- PUBLISHED: 6 00/422 
US - CL - CURRENT : 600 / 422 

REPRESENTATIVE-FIGURES: 8A 
ABSTRACT : 

The present invention is a technique of, and system for, imaging vascular anatomy 
over distance considerably greater than the maximum practical field of view of a 
magnetic resonance imaging system while using substantially one contrast agent 
injection. The technique and system of the present invention acquires image data of 
a plurality of image volumes which are representative of different portions of the 
patient's body. The image data of each image volume includes i mage data which is 
representative of the center of k- space . The acquisition of image data which is 
representative of the center of k-space is correlated with a concentration of 
contrast agent in the artery (ies) residing in the image volume being substantially 
greater than the concentration of contrast agent in veins and background tissue 
adjacent to the artery (ies) . This provides preferential enhancement of arteries 
relative to adjacent veins and background tissue for each acquisition, wherein each 
acquisition is representative of a different portion of the arterial system (e.g., 
abdominal aorta, femoral, popliteal, and tibial arteries) . 
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ABSTRACT : 

MRA data is acquired from a large region of interest by translating the patient 
through the bore of the MRI system as a three-dimensional MRA data set; are acquired. 
Patient table movement is controlled to track a bolus of contrast agent as it passes 
through the region of interest . Fluoroscopic images may be acquired during the scan 
to enable accurate bolus tracking. A seamless image of the entire region of interest 
is reconstructed. 

CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] This application is based on U.S. Provisional Patent Application Serial No. 
60/282,555 filed on Apr. 9, 2001 and entitled "Method For Acquiring MRI Data From A 
Large Field Of View Using Continuous Table Motion." 
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TITLE: Sampling and reconstruction of signals and images including MR images of 
multiple regions 



Abstract Text (1) : 

A method and apparatus are provided for acquiring and reconstructing an image . The 
method includes the steps of obtaining prior knowledge of the image/ possibly by 
coarse sampling of the image and using the obtained prior knowledge of the image to 
identify relative locations of structures having relatively high contrast edges . The 
method further includes the steps of prescribing a set of k-space locations based 
upon the relative locations of the structures in order to achieve comparable 
eigenvalues of a reconstruction matrix and sampling the k-space at the prescribed 
k-space locations to obtain k-space sample data. The k-space sample data are 
decomposed into background data and edge data. The background data are Fourier 
transformed to reconstruct a background image component. Similarly, subsets of the 
edge data are Fourier transformed and the reconstruction matrix is used to form a 
linear combination of these Fourier transformations in order to reconstruct an edge 
image component. Finally, the background image component and the edge image 
component are combined to generate a final image . 

Brief Summary Text (2) : 

The field of the invention relates to sampling and reconstruction of signals and 
images including MR images of multiple regions. 

Brief Summary Text (4) : 

MICROFICHE APPENDIX I contains one page of microfiche with 51 frames and is a 
printout of source code for executing steps of an illustrated embodiment of the 
invention" . 

Brief Summary Text (5) : 

MR imaging is most commonly performed with 2D FT or 3D FT techniques, which require 
only a modest amount of prior knowledge about the object: 1) the object is assumed 
to be completely contained within a finite field of view ( FOV ) ; 2) the image is 
assumed to be band-limited in frequency space ; i.e. its significant power spectrum 
does not extend beyond some maximum spatial frequency . If these two assumptions are 
valid, the usual sampling theorem (due to Whittaker, Kotel'nikov, and Shannon, or 
WKS, (1-4)) states that the image can be reconstructed by an inverse FT of a finite 
number of discrete samples of the image ' s k- space representation. The spacing 
between the sampled k- space points is determined by the first assumption (the 
dimensions of the FOV ) , and the limits of the k -space sampling pattern are 
determined by the second assumption (the maximum spatial frequencies in the image ) . 
This method is widely used because: 1) it is based on relatively weak assumptions 
about the image , 2) it is relatively easy to acquire the Fourier -encoded signals 
stipulated by the WKS theorem, and 3) H5?.H®. reconstruction can be performed 
efficiently with a fast FT (FFT) . 

Brief Summary Text (6) : 

In contrast-enhanced carotid artery „ i^SiQ&t^ interventional imaging, functional 
:y?.?SA?.3„'.„ cardiac iTOQing_,_ and a number of other applications, "the utility of MRI is 
limited by the speed with which the l^.-s^ace^ data can be measured. Two general 
strategies have been used to shorten A^aae acquisition time. 1) 5£^lient pulses with 
shorter rise times and/or larger amplitudes have been used in order to shorten the 
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time required to gather a complete WKS data set . Unfortunately, gradient: ramp rates 
and strengths are now approaching values at which neuromuscular stimulation can 
compromise patient safety and comfort. 2) More stringent assumptions can be made 
about the image in order to reduce the number of signals necessary to reconstruct 
it. Some of these "constrained imaging 11 methods simply apply WKS sampling with 
stronger assumptions in order to increase the spacing between k- space points 
(reducing the FOV ) or to reduce the k-space sampling limits (reducing image 
resolution) . More novel approaches have utilized prior knowledge to express the 
image as a superposition of a small number of non-Fourier basis functions. The 
image ' s projections onto these basis functions are computed from a reduced set of 
Fourier -encoded signals, or they are measured directly by performing non-Fourier 
encoding . 



Brief Summary Text (7) : 

In this invention, we take a different approach. We generalize the WKS sampling 
theorem so that it can be applied to images which are supported on multiple regions 
within the FOV . By using this "multiple region MR" (mrMR) sampling theorem, such 
images can be reconstructed from a fraction of the k-space samples required by the 
WKS theorem. Image reconstruction is performed with FFTs and without any noise 
amplification, just as in conventional FT MRI . In addition, we show how the method 
can be applied to a broader class of images having only their high contrast edges 
confined to known regions of the FOV . If this kind of prior knowledge is available, 
k-space can be sampled sparsely, and scan time can be reduced. The next section 
describes the theoretical framework of the mrMR approach. Then, the method is 
illustrated with simulated data and with experimental data from a phantom. Finally, 
we describe how the method was used to reduce the time of first-pass Gd-enhanced 3D 
carotid MRA so that it could be performed without bolus timing. 

Brief Summary Text (9) : 

A method and apparatus are provided for acquiring and reconstructing an image . The 
method includes the steps of obtaining prior knowledge of the image , possibly by 
coarse sampling of the image and using the obtained prior knowledge of the image to 
identify relative locations of structures having relatively high contrast edges . The 
method further includes the steps of prescribing a set of k-space locations based 
upon the relative locations of the structures in order to achieve comparable 
eigenvalues of a reconstruction matrix and sampling the k-space at the prescribed 
k-space locations to obtain k- space sample data. The k- space sample data are 
decomposed into background data and edge data . The background data are Fourier 
transformed to reconstruct a background image component. Similarly, subsets of the 
edge data are Fourier transformed and the reconstruction matrix is used to form a 
linear combination of these Fourier transformations in order to reconstruct an edge 
image component. Finally, the background image component and the edge image 
component are combined to generate a final image . 

Brief Summary Text (10) : 

Traditional Fourier MR imaging utilizes the Whittaker-Kotel ' nikov-Shannon (WKS) 
sampling theorem. This specifies the spatial frequency components which need be 
measured to reconstruct an image with a known field of view ( FOV) and band- limited 
spatial -frequency contents. In this paper, we generalize this result in order to 
find the optimal k-space sampling for images that vanish except in multiple, 
possibly non-adjacent regions within the FOV . This provides the basis for "multiple 
region MRI" (mrMRI), a method of producing such images from a fraction of the 
k -space samples required by the WKS theorem. Image reconstruction does not suffer 
from noise amplification and can be performed rapidly with fast Fourier transforms, 
just as in conventional FT MRI . The mrMRI method can also be used to reconstruct 
images that have low spatial -frequency components throughout the entire FOV and high 
spatial frequencies (i.e. edges ) confined to multiple small regions. The greater 
efficiency of mrMR sampling can be parlayed into increased temporal or spatial 
resolution whenever the imaged objects have signal or " edge" intensity confined to 
multiple small portions of the FOV. Possible areas of application include MR 
angiography (MRA), interventional MRI, functional MRI, and spectroscopic MRI. The 
technique is demonstrated by using it to acquire Gd-enhanced first-pass 3D MRA 
images of the carotid arteries without the use of bolus-timing techniques. 

Drawing Description Text (2) : 
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FIG. 1 depicts a schematic outline of an illustrated image reconstruction method in 
accordance with an embodiment of the invention; 

Drawing Description Text (4) : 

FIG. 3 depicts an image supported on three cells that may be processed using the 
method of FIG. 1; 

Drawing Description Text (5) : 

FIG. 4 depicts the application of the method of FIG. 1 and FIG. 2 to images of a 
phantom; 

Drawing Description Text (6) : 

FIG. 5 depicts further details of the application of the method to the phantom 
images ; 

Drawing Description Text (7) : 

FIG. 6 depicts an image of a volunteer processed by the method of FIG. 1 and FIG. 2; 



Drawing Description Text (8) : 

FIG. 7 depicts additional detail of the image of FIG. 6; 
Drawing Description Text (10) : 

FIG. 9 is a block diagram of apparatus that may be used to practice the method of 
FIG. 1 and FIG. 2. 

Detailed Description Text (3) : 

The mrMRI method of an illustrated embodiment of the invention is schematically 
outlined in FIGS. la-e. In FIG. la, the image to be reconstructed (top panel) 
vanishes except in three cells that covered the FQV (bottom panel) . In FIG. lb, the 
composite sampling pattern used to reconstruct the mrMR image from a small number of 
Fourier components is shown. This pattern is composed of a number of sparse grids, 
each offset by a different amount from the center of k- space . In this example, the 
data sampled at all of these points comprise just 1/12 of the WKS data set needed 
for reconstruction of the entire FQV with the same resolution. In FIG. lc is shown 
the "reduced FQV" images, S.sub.a (x,y), which were produced by inverse FT of data 
on each of the sparse grids in panel b (FIG. lb) . Each of the images is a linear 
combinations of all the cell images, I.sub.j (x,y). Only the magnitudes of these 
complex images are shown. In FIG. Id, the contents of each cell, found by linearly 
combining the "reduced FQV" images in panel c is shown. In FIG. le, the mrMR image, 
created by tiling the FQV with the cell images in panel d is shown. This image is 
windowed identically to the "exact" image in panel a. 

Detailed Description Text (4) : 

A flow chart of the process depicted in FIGS, la-e may be as shown in FIG. 8. 
Apparatus 100 that may be used to practice the process shown in FIGS. 1 and 8 may be 
as shown in FIG. 9. Reference shall be made to these figures as appropriate to 
understanding the invention. 

Detailed Description Text (5) : 

Consider an arbitrary image in the xy-plane (top panel of FIG. la) . Cover this plane 
with an infinite array of rectangular cells with any convenient dimensions (. DELTA. x 
and . DELTA. y) along the x and y axes, respectively (bottom panel of FIG. la) . If the 
image were supported on only one cell, the WKS sampling theorem would dictate that 
k -space be sampled by the data acquisition device 102 of FIG. 9 on a single coarse 
grid with spacing equal to . DELTA . k . sub . x =2 .pi . / . DELTA . x and . DELTA. k . sub . y 
=2 .pi . / .DELTA. y along the k.sub.x and k.sub.y axes, respectively. It can be shown 
that an image supported on C such cells can be reconstructed from sparse data on C 
(or more) such sparse grids, each grid being offset by a different amount from the 
center of k- space (FIG. lb) . 

Detailed Description Text (6) : 

First, let I(x,y) be the image to be reconstructed, and denote its Fourier transform 
by I (k . sub .x, k . sub . y) : ##EQU1## where all integrations extend from minus infinity to 
plus infinity, unless specified otherwise. Define the "comb" sampling function, 
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III(k), to be ##EQU2## Next, define the function S.sub.a (k . sub .x, k . sub . y) to be 
##EQU3## where k.sub.a . ident . (k . sub . ax, k . sub . ay) is the offset between the a.sup.th 
sampling grid and the origin of k- space ; a=l, . . . , A; and A is the number of 
coarse grids that are superposed to form a composite sparse sampling pattern. 
Equations (2) and (3) show that S.sub.a (k . sub .x, k . sub . y) simply represents the 
image ' & Fourier components on the a.sup.th grid and can be determined experimentally 
by measuring MR signals that are Fourier -encoded with the appropriate wave vectors. 
Taking the inverse FT of Eq. (3) produces (FIG. lc) ##EQU4## Let the integers (m, n) 
label each cell in the xy-plane with respect to any convenient cell that is labeled 
{0, 0) and is centered at (x.sub.O, y.sub.O) (FIG. la). Let us define I . sub . j (x,y) 
to be the function produced by translating the contents of the j . sup . th cell, 
located at cell coordinates (m.sub.j, n.sub.j), to the cell at (0, 0): 

Detailed Description Text (10) : 

and ##EQU8## In other words, each experimentally-determined function S.sub.a (x, y) 
is a linear combination of all the cell images I.sub.j (x,y) . Notice that the 
contribution of each cell is weighted by a constant complex number, which depends on 
the offset of the coarse grid from the center of k-space (k.sub.a) and on the 
locations of the cells (m.sub.j and n.sub.j) 

Detailed Description Text (11) : 

Up to this point, no assumptions have been made about the image . Now let us apply a 
more flexible version of the WKS constraint on the FOV . Assume that only C of the 
cell images are non-zero; i.e. all of the remaining elements of I are zero. In this 
case the A. times . . inf in . matrix M can be replaced by an A. times. C matrix, and the 
. inf in . -dimensional vector I can be replaced with a C-dimensional vector, with j 
varying over only the non-vanishing cells. If we have measured data on A coarse 
grids (A . gtoreq . C) , we are left with a system of A equations in C unknowns. The C 
non-vanishing functions I.sub.j can be found by inverting the A linear equations of 
the form of Eqs . (6) and (7). The solution is ##EQU9## j . epsilon . non-zero cells, (9) 
or in matrix notation 

Detailed Description Text (13) : 

Notice that we have only imposed a spatial constraint on the image ' s region of 
support. We have not made any assumptions about the image 1 s spatial frequency 
content within the signal -containing cells. This latter assumption determines the 
k -space limits (and therefore the total size) of each of the A sparse grids, and 
ultimately determines the resolution of the image within the cells. 

Detailed Description Text (14) : 
B - Optimal Sampling 

Detailed Description Text (15) : 

Equations (9) and (10) can be used to reconstruct the image from experimental data 
on A coarse k- space grids with appropriate offsets. Notice that the matrix M depends 
on the values of these offsets, as well as on the locations of the supporting cells 
(Eq. (8)) . In the absence of noise, any offsets leading to a non-singular M can be 
used to accurately reconstruct the image, as long as the number of coarse sampling 
grids is at least as great as the number of non-zero cells in the image (i.e. 
A. gtoreq. C) . However, in the presence of noise, the conditioning of M (or 
equivalently M.sup.+ M) determines the amount of noise amplification during 
reconstruction. Specifically, noise in the experimental data will be amplified 
during the reconstruction process to the extent that M is ill-conditioned (i.e. to 
the extent that some of the eigenvalues of M.sup.-i- M are much smaller than others) . 
Therefore, one expects that image noise will be minimized if the k-space offsets are 
chosen so that these eigenvalues are all equal . 

Detailed Description Text (16) : 

To prove this, we use Eqs. (9) and (10) to calculate the mean squared noise 
propagating from the experimental data into the reconstructed image I (x, y) . Assume 
that the root -mean -squared noise in I (k . sub .x, k . sub . y) is .epsilon. and that there 
is no correlation between the noise at different k-space locations. Then, it follows 
from Eqs. (9) and (10) that the mean squared noise in the 2-^9^ is: ##EQU10## where 
the brackets denote the statistical average, the integral extends over the C cells 
that support I, and each coarse k-space grid has the dimensions N.sub.x 
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. times. N. sub. y. The trace on the right side is equal to ##EQU11## where .mu..sub.i 
denotes the i.sup.th eigenvalue of the C.times.C Hermitian matrix M.sup.+ M. Notice 
that Eq. (8) implies that these eigenvalues must lie on a hyperplane in .mu. -space 
defined by ##EQU12## Straight -forward application of multivariable calculus shows 
that the trace in Eq. (12) has a global minimum on this hyperplane when all 
.mu..sub.i =A (i=l, . . . , C) . Therefore, as we expected, the level of noise in an 
mrMR image is minimized when the A coarse sampling pattern offsets are chosen such 
that these eigenvalues are all equal. Notice that for this optimal sampling pattern, 
##EQU13## where NEX . ident . A/C is the ratio of the number of acquired k- space points 
to the minimal number. 

Detailed Description Text (17) : 

Notice that the eigenvalues of M.sup.+ M are equal to one another if and only if 
this matrix is equal to A times a C.times.C identity matrix. Using the definition of 
M (Eq. (8)), M.sup.+ M has this form if and only if ##EQU14## for all unequal j and 
j', where d.sub.jj, is the displacement vector between the centers of any two 
supporting cells, located at (m.sub.j, n.sub.j.) and (m.sub.j, n.sub.j.). Equation 

(14) shows explicitly how the optimal offsets of the sampling pattern (k.sub.a) are 
determined by our prior knowledge of the locations of the supporting cells 

(d.sub.jj 1 ). Notice that Eq. (14) only depends on the relative locations of the 
supporting cells. Therefore, the optimal sampling pattern is not changed if the cell 
array is translated across the FOV in order to most efficiently contain the 
intensity distribution of a particular image . 

Detailed Description Text (18) : 

To summarize, we have shown that the offsets of the A coarse k- space grids should be 
chosen so that Eq. (14) is satisfied for all pairs of supporting cells. If this can 
be achieved, image reconstruction (Eqs. (9) and (10)) will be perfectly 
well -conditioned, and image noise will be minimized. In fact, the level of noise 
will be the same as in a conventional (WKS sampled) FT image having the same 
resolution over an equal area of support (i.e. over a rectangular FOV with area 
equal to C. DELTA. x . DELTA. y) . This last statement is a consequence of the following 
two facts: 1) the level of noise in an mrMR image does not depend on the positions 
of the supporting cells; 2) if the supporting cells are rearranged into a 
rectangular array (e.g. into a single row of adjacent cells), the WKS sampling 
pattern is an optimal (mrMR) sampling pattern. The first fact follows immediately 
from Eq. (11) and the fact that Tr[(M.sup.+ M).sup.-1 ] is independent of the 
supporting cells' positions; namely, for an optimal sampling pattern, it is equal to 
C/A or 1/NEX. To understand the second statement, imagine that the supporting cells 
form an N.sub.x . times .N. sub . y rectangular array. The WKS sampling pattern for this 
"reduced FOV " can be broken down into a composite of A=N.sub.x N.sub.y coarse 
k- space grids, having the offsets: ##EQU15## where p. sub. a =0, 1, ... , N.sub.x -1 
and q.sub.a =0, 1, ... , N.sub.y -1. For these offsets, the mrMR matrix element 
M.sub.aj in Eq. (8) is e . sup . -2 . pi . ip . sbsp . a . sup . m. sbsp . j . sup . /N . sbsp .x 
e . sup . -2 .pi . iq . sbsp . a . sup .n . sbsp . j . sup . /N. sbsp .y . Each of these factors is a 
component of a discrete Fourier function. It follows from their completeness and 
orthogonality that M.sup.+ M=N.sub.x N.sub.y I=A I, where I is the identity matrix. 
Therefore, the WKS sampling pattern leads to an M.sup.+ M matrix having all 
eigenvalues equal to A, proving the assertion that it is also an optimal mrMR 
sampling pattern. Therefore, we have proved that the noise level in an mrMR image is 
the same as that for a WKS image of a rectangular region formed by the rearrangement 
of the supporting cells. 

Detailed Description Text (20) : 

We now show how the system 100 may use a variant of mrMRI to reconstruct an image 
that has non-zero intensity throughout the entire FOV but has high contrast edges 
confined to multiple regions within the FOV. Suitable filters can be used to 
decompose such an image, into the sum of a slowly undulating "background" component 
and an "edge" component that is only supported on small portions of the FOV- The 
former component can be reconstructed from data in a small "keyhole" at the center 
of k- space, and the latter component can be reconstructed from sparse .„]?." space data 
in an optimal mrMR sampling pattern. FIG. 2 outlines this filtered mrMRI technique. 

Detailed Description Text (21) : 

FIG. 2a, shows the image to be reconstructed (left panel) consisted of a slowly 
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undulating background component, superposed on an edge component that was supported 
on only three cells (right panel) . FIG. 2b, shows the k- space representation of the 
Hamming filter which may be used to decompose the image into background and edge 
components. FIG. 2c. shows the background component (right panel) and its FT (left 
panel), which may be used to define the product of the FT of the exact image (panel 
a) and the Hamming filter (panel b) . FIG. 2d shows the edge component (right panel) 
and its FT (left panel), which may be used to define the product of the FT of the 
exact image (panel a) and the complement of the Hamming filter (panel b) . FIG. 2e 
shows the background component (right panel) that may be reconstructed by densely 
sampling the center of the image's k- space using a keyhole sampling pattern (left 
panel), applying the low-pass filter (panel b) to these data, and taking the inverse 
FT of this filtered data. FIG. 2f shows the edge component (right panel) that may be 
reconstructed by sampling the image's k- space on each of the sparse grids forming an 
mrMR sampling pattern (left panel) , applying the complementary high-pass filter to 
these data, and then applying the mrMR technique (FIG. 1) . FIG. 2g shows the mrMR 
image, created by adding the reconstructions of the background and edge components 
(panels e and f ) . This image was reconstructed from just 1/9 of the WKS k-space data 
required to create an image of the entire FOV with the same resolution. This image 
is windowed identically to the exact image in panel a. FIG. 2h shows the error image 
formed by taking the magnitude of the difference between the mrMR image (panel g) 
and the exact (WKS) image (panel a) . This image has been windowed to accentuate the 
faint discrepancies. 

Detailed Description Text (22) : 

First, note that the FT of any image can be decomposed into complementary components 
(FIGS. 2c and 2d) that are low-pass filtered and high-pass filtered, respectively: 

Detailed Description Text (23) : 

where I (k . sub .x, k . sub . y) is the FT of the image and F (k . sub .x, k . sub . y) is the 
low-pass filter (FIG. 2b). If F (k . sub . x, k . sub . y) is chosen wisely, B(x,y) will 
resemble I(x,y) except in the vicinity of high contrast edges which will be modified 
by blurring and truncation artifact (FIG. 2c), and E(x,y) will vanish everywhere 
except in the vicinity of the edges of I(x,y) (FIG. 2d) . By definition, the 
background component, B(x,y), can be reconstructed from data in a k- space "keyhole" 
that is bounded by the cut-off of F (k . sub . x, k . sub . y) (FIG. 2e) . Specifically, data 
can be acquired on the dense WKS grid in this region, multiplied by 
F (k . sub .x, k . sub.y) and then subjected to an FFT. The result is B(x,y), which 
comprises a low resolution version of the filtered image . The mrMR method can be 
used to reconstruct the edge component, E (x, y) , from E (k . sub .x, k. sub . y) measured at 
optimal sparse k- space locations (FIG. 2f ) . The required values of 
E (k . sub .x, k . sub . y) can be determined experimentally by measuring the corresponding 
components of I (k . sub . x, k . sub . y) and multiplying them by the high pass filter 
1-F (k . sub .x, k . sub . y) . Finally, the image I(x,y) is found by summing the 
reconstructed background and edge components: I (x, y) =B (x, y) +E (x, y) , as shown in FIG. 
2g. 

Detailed Description Text (24) : 

Neither B(x,y) nor E (x, y) will be degraded by noise amplification because both are 
reconstructed from optimal k- space sampling patterns, determined by the WKS and mrMR 
sampling theorems, respectively. E(x,y) could be degraded by aliasing if the cut-off 
frequency of F(k.sub.x, k. sub.y) is too low, causing ringing and blurring patterns 
in E(x,y) to extend so far from the edges in I(x,y) that they reach beyond the 
boundaries of the cells used in the mrMR reconstruction. However, this problem can 
be ameliorated by: 1) using a suitably high cut-off frequency in F (k . sub .x, k . sub . y) ; 
2) designing the shape of F(k.sub.x, k. sub.y) so that the truncation artifact is 
largely confined to the vicinity of the edges (e.g. a Hamming filter); 3) using 
cells that are sufficiently large to contain most of the remaining truncation 

Detailed Description Text (25) : 

In the special case of an ..A™?.^ confined to a single region, the filtered 

mrMR method may be similar to that of certain prior art techniques. Such -"i^S.? 8 ma Y 
be reconstructed from k-spac_e samples on a single sparse WKS grid, superposed on a 
pattern of dense sampling in the center of • However, the prior art 

techniques used a step_-l.ike_ function for F(k.sub.x, k. sub.y). This choice of filter 



6 of 20 



10/27/2003 9:01 AM 



Record Display Form http://westbrs:8002/biiV gate. exe?f=dn^s. . .c_2=&p_doc_3=&p_doc_4=&p_doc_5=&p_doc_6= 



was suboptimal because F (x, y) in this case was a sine function, and edges in I (x, y) 
produced non-localized ringing in E(x,y). This tended to increase the region of 
support of E(x,y) and reduced the efficiency of the method. 

Detailed Description Text (27) : 

Suppose that an image is known to be supported on only three cells (FIG. 3) and that 
tne ^ m age is to b e reconstructed from data on three coarse grids. FIG. 3 shows image 
that is supported on only three cells, which are separated from each other by 
displacement vectors d.sub.ij. In this case, Eq. (14) is equivalent to the following 
three conditions: 

Detailed Description Text (29) : 

First, consider the special case in which the three cells happen to coalesce into a 
single rectangular region parallel to the x -axis . Then, we have d. sub. 12 =d.sub.23 
=(. DELTA. x, 0) and d. sub. 31 = ( -2 . DELTA . x, 0), and algebra shows that all solutions 
to Eq. (16) must satisfy: k.sub.lx -k.sub.3x =1/3 . DELTA. k . sub . x and k.sub.2x 
-k.sub.3x =2/3 .DELTA. k. sub. x. The three superposed coarse grids with these offsets 
comprise the usual WKS sampling pattern for a single region with dimensions 
3 .DELTA. x. times . .DELTA. y. In other words, the general condition for optimal mrMR 
sampling (Eq. (16)) reduces to the usual sampling theorem in the special case when 
the three cells coalesce into a single rectangular region. This illustrates the 
general statement which was proven at the end of IIB: WKS sampling patterns belong 
to the larger set of optimal mrMR sampling patterns. 

Detailed Description Text (30) : 

Now, suppose that the three cells are at non-collinear locations (FIG. 3) . With the 
help of the fact that d. sub. 12 +d.sub.23 +d.sub.31 =0, it can be shown that the only 
solutions of Eq. (16) are: ##EQU16## where n.sub.i (i=l, . . . , 4) are any integers 
and D is the matrix: ##EQU17## For example, consider the image in FIG. la, which is 
supported on only three cells of a 6. times. 6 cell array. In this case, d. sub. 12 
= (4 . DELTA. x, 0) and d. sub. 23 = (0, 4 . DELTA . y) , and the solutions in Eq. (17) are 
k.sub.l -k.sub.3 = (1/3 .DELTA. k. sub. x, 5/6 .DELTA. k. sub .y) and k.sub.2 -k.sub.3 
= (1/6 . DELTA. k. sub .x, 2/3 . DELTA. k . sub . y) when we choose n.sub.i =1, n.sub.2 =3, 
n.sub.3 =-1, and n.sub.4 =-3. FIG. lb shows the optimal composite sampling pattern, 
created by superposing these three coarse grids and taking k.sub.3 =0 for 
convenience. Other optimal sampling patterns can be created by substituting 
different values for n.sub.i. For example, the choices of n.sub.i =n.sub.2 =1 and 
n.sub.3 =n.sub.4 =-3 give k.sub.l -k.sub.3 =( 1/3 . DELTA . k . sub . x, 1/3 . DELTA . k . sub . y) 
and k.sub.2 -k.sub.3 = (2/3 .DELTA. k. sub. x, 2/3 . DELTA . k . sub . y) 

Detailed Description Text (31) : 

The mrMR method may be tested by using it to reconstruct the 192 . times . 192 image in 
FIG. la from simulated data in an optimal mrMR sampling pattern. Simulated data may 
be generated by Fourier transformation of the image in FIG. la and then sampled 
sparsely according to the sampling pattern in FIG. lb. The number of sampled data 
was 1/12 of the data needed for a conventional FT reconstruction of the entire FOV 
with the same resolution. These simulated data may be substituted into the right 
side of Eq. (9) or (10) in order to reconstruct the mrMR version of the image (FIG. 
le) . 

Detailed Description Text (32) : 

The filtered mrMR method may be tested by using it to reconstruct the 192 . times . 192 
image in FIG. 2a from simulated data. The filter, F(k.sub.x, k.sub.y), may be chosen 
to be a Hamming filter that was cut off along the edges of a 32. times. 32 square in 
the center of k - space (FIG. 2b) . The effect of this filter is shown in FIGS. 2c and 
2d, which depict the background and edge components of the image in FIG. 2a. By 
definition, the background component is band- limited to a 32. times. 32 region of 
k- space . Therefore, it can be reconstructed exactly from a 32. times. 32 array of data 
on "the" WKS grid (FIG. 2e) . Notice that the intensity in the edge image in FIG. 2d is 
confined to the same three cells that support FIG. la. Therefore, this .^9^..A!?3.3 e „ 
should be accurately reconstructed from the sparse mrMR sampling patterns derived 
for FIG. la (namely, the sampling pattern in FIG. lb or FIG. 2f ) . The composite 
sampling pattern may be created by combining the "keyhole" of data required to 
reconstruct the background component (FIG. 2e) with the ®^3®. component sampling 
scheme (FIG. 2f ) . These combined data comprised only 1/9 of the complete WKS set. 
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Simulated data may be generated by Fourier transforming the image in FIG. 2a and 
then sampled sparsely at the locations in FIGS. 2e and 2f . The reconstructed 
background and edge images (FIGS. 2e and 2f) may be added in order to generate the 
mrMR version of the image (FIG. 2g) . Then, an error image (FIG. 2h) may be formed by 
taking the magnitude of the difference between the exact image (FIG. 2a) and the 
mrMR image (FIG. 2g) . 

Detailed Description Text (34) : 

To demonstrate the validity of the technique described above, the mrMRI technique 
was used to image a carotid artery phantom (FIGS. 4 and 5) consisting of two 
branching plastic tubes having an inner diameter equal to 5 mm. FIG. 4a shows a low 
resolution "scout" image of the carotid artery phantom, with a superposed 4. times. 8 
array of rectangular cells. FIG. 4b shows the same scout image after the cell array 
was translated and "sheared" in the y -direction so that all the vessels fell within 
just two cells. Note that translation and shearing of the cell array (possibly on a 
slice-by slice basis) does not change the optimal sampling pattern dictated by mrMR. 



Detailed Description Text (35) : 

FIG. 5a shows the "exact" image of a typical axial slice of the phantom, generated 
by FT of 16,384 ( 128 . times . 128 ) phase -encoded signals. FIG. 5b shows the k- space 
locations of Fourier components that are optimal for mrMR reconstruction of all 
slices in the phantom. These comprise less than 1/8 of the WKS data set used to 
produce the exact image in panel a. FIG. 5c shows the mrMR image of the axial slice 
in panel a, reconstructed from the signals at the k- space locations in panel b. This 
image is windowed identically to panel a. FIG. 5d shows the "error" image produced 
by subtracting the exact image in panel a from the mrMR reconstruction in panel c 
and windowing the result to increase the conspicuity of the faint discrepancies. 
FIG. 5e shows MIPs created from all of the exact axial images . FIG. 5f shows MIPs 
created from the mrMR reconstructions. This image is windowed identically to panel 
e. FIG. 5g shows an image of the errors in the mrMR MIPs, produced by subtracting 
the exact MIPs in panel e from the mrMR MIPs in panel f and windowing the result to 
accentuate the faint discrepancies. 

Detailed Description Text (36) : 

The tubes were filled with a concentrated gadolinium solution (gadodiamide 
( Omni scan ) , Nycomed, Inc., Wayne, Pa.) and were immersed in a bath of very dilute 
gadolinium solution. The phantom was placed in the quadrature bird-cage head coil of 
a 1.5 T whole body scanner (Signa, General Electric Medical Systems, Milwaukee, 
Wis.). A 3D spoiled gradient echo pulse sequence (17) was used with TR=4 . 1 ms, 
TE=1.6 ms, flip angle=20. degree. . This pulse sequence selectively excited a coronal 
slab that had RL. times .AP . times .SI dimensions of 160 . times . 160 . times . 160 and 
contained the "arteries". The slab was phase -encoded in the transverse plane and 
frequency - encoded along the SI direction to produce a 128 . times . 128 . times . 128 data 
array, resulting in an isotropic 1.25 mm resolution. These 16,384 k- space signals 
were obtained in 67 seconds. 

Detailed Description Text (37) : 

Next, a small subset of these signals was used to reconstruct mrMR images of each 
axial section. The filtered mrMR technique was applied in order to account for the 
bath's signal intensity, which varied slowly across the entire FOV . The background 
component of each image was defined by Eq. (15), with F (k . sub .x, k . sub . y) equal to a 
Hamming filter with a cut-off at the edges of a 32. times. 32 square in the center of 
k- space . This background image was reconstructed by conventional FFTI of the signals 
that densely sampled this region, after they had been multiplied by the Hamming 
filter. The complementary edge image was expected to vanish except in the vicinity 
of the vessels, and was reconstructed from k -space samples in the sparse mrMR 
pattern dictated by the locations of the vessel-containing cells. In an actual 
clinical situation, a scout scan of the neck would have been used to ascertain the 
relative positions of the vessels, as demonstrated in the next section. In the case 
of the phantom, the vessels were located on low resolution (32 . times . 32 ) axial 
sections of the above -described coronal slab. The FOV of these axial sections was 
covered by a 4. times. 8 array of rectangular cells (FIG. 4a) with dimensions: 
. DELTA. x=40 mm and . DELTA . y=2 0 mm. This array was translated and sheared so that the 
vessels were contained in just two cells (i.e. C=2) , which were separated by one 



8 of 20 



10/27/2003 9:01 AM 



Record Display Form * http://westbrs: 8002/bin/gate. exe?f=do^&s. . . c_2=&p_doc_3 =&p_doc_4=&p_doc_5 =&p_doc_6= 



empty cell along the x -direction (FIGS. 4a and 4b) . As shown below, the contents of 
the sheared cells can be reconstructed without noise amplification from the optimal 
mrMR sampling pattern for unsheared cells. An mrMR edge image for each axial section 
was reconstructed from a sparse subset: of the acquired signals, which included the 
k- space points on two coarse grids, i.e. A-2 . The mrMR sampling criterion in Eq. 
(14) dictates that the offsets of these grids (k.sub.l and k.sub.2 k) satisfy: 

Detailed Description Text (39) : 

where J is an odd integer. Because the displacement between the centers of the 
unsheared edge -containing cells was d= (2 . DELTA . x, 0), Eq. (19) implies that k.sub.lx 
-k.sub.2x =J/4 .DELTA. k. sub. x, where . DELTA . k . sub . x =2 . pi . / . DELTA . x and J is any odd 
integer. This constraint was satisfied with J=l, k.sub.l =( 1/4 . DELTA. k . sub . x, 0), 
and k.sub.2 =0, which corresponds to the two coarse grids in FIG. 5b. 

Detailed Description Text (40) : 

Parenthetically, it is interesting to notice what would have happened if the cells 
had been adjacent (d= ( .DELTA. x, 0) ) . In that case, Eq. (19) would have dictated that 
k.sub.lx -k.sub.2x =J/2 .DELTA. k. sub. x, where J is again an odd integer. For all 
values of J, this implies that the two coarse grids must be offset by 
1/2 .DELTA. k. sub. x with respect to one another, thereby forming the usual WKS 
sampling pattern for a single rectangular region with dimensions 

2 .DELTA. x. times .. DELTA. y. This example is another illustration of the fact which was 
proven in general at the end of I IB: the optimal mrMR sampling patterns include to 
the usual WKS sampling when the cells coalesce into a single rectangular region. 

Detailed Description Text (41) : 

In order to reconstruct the edge component of the image, the data on each of the 
32. times. 16 coarse grids was multiplied by the complement of the Hamming filter, and 
each of these filtered data sets was inverse Fourier transformed. The resulting 
functions S . sub . a (x,y) were used to compute images of the two sheared 
edge -containing cells with 1.25 mm. times . 1 . 25 mm resolution. The resulting edge 
image was added to the previously-described background image to create the mrMR 
image (FIG. 5c) . This should be compared to the exact image (FIG. 5a) , which was 
reconstructed by FT of the 16,3 84 data in the complete (128 . times . 128) WKS sampling 
pattern. Notice that the mrMR image was reconstructed from the 1,984 data in the 
composite sampling pattern (FIG. 5b), which comprised less than 1/8 of the complete 
data set . FIG. 5d shows the error image that was equal to the magnitude of the 
difference between the exact image (FIG. 5a) and the mrMR image (FIG. 5c) . Finally, 
maximum intensity projection (MIP) images were created from the exact and mrMR 
reconstructions (FIGS. 5e and f ) , and differences between these MIPs were displayed 
in corresponding error images (FIG. 5g) . 

Detailed Description Text (42) : 

As mentioned above, for some axial images, the cell array had to be sheared along 
the y - axi s so that the vessels were contained in just two cells (FIG. 4b) . However, 
it is easy to see that the contents of the sheared cells can still be reconstructed 
from the mrMR sampling pattern derived for the unsheared cells. In fact, this can be 
done whenever each non-zero cell occupies a different cell column and the k.sub.a 
offsets all lie on the k.sub.x axis . In this situation, the sparse sample spacing in 
the k.sub.y direction is determined only by the y-extent of the image intensity in 
each pixel column, independent of the cell width. Therefore, the 2D mrMR 
reconstruction is equivalent to a ID mrMR reconstruction in the x-dimension and a 
conventional ID "reduced FOV " reconstruction in the y-dimension. As long as the 
image intensity in each column does not extend over a segment greater than the 
height of a single cell (. DELTA. y), it can be positioned in a continuous manner 
anywhere within that particular column, and the resulting cell array will appear to 
be sheared in the y-dimension. The additional flexibility of this shearing strategy 
allows one to customize the cell shapes on a slice by slice basis without affecting 
the °P5i?i?;L mr MR sampling pattern. This reduces both the height of the cells needed 
to contain the i^S.5r..!..?....?^9?.?... anc * tne number of ^--sgace^ points needed to successfully 
reconstruct the r mage :. it should be noted that this shearing method can also be 
applied to certain arrangements of supporting cells in which some cell columns 
contain more than one supporting cell; however, such cases were not encountered in 
the experiments described here. 
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Detailed Description Text (44) : 

The mrMR method was used to perform first -pass Gd- enhanced 3D MRA of the carotid 
arteries of a patient volunteer. Informed consent was obtained from the patient, a 
68 year-old male who was scheduled to undergo an infused brain MR I . After the 
uninfused portion of the brain exam was completed f the patient underwent an 
unenhanced 2D TOF neck MRA with the same neurovascular coil used for the brain 
study. These images were obtained for comparison to the subsequent Gd-enhanced mrMRA 
images . 2D TOF was performed with a spoiled gradient; echo pulse sequence with TR=23 
ms, TE=4 . 9 ms, flip angle=60 . degree . , first order flow compensation, and fractional 
echo. Each axial section was 1 . 5 mm thick and was in-aged with a 256 . times . 224 matrix 
across a 180 mm. times. 18 0 mm FOV . Thus, the overall spatial resolution was 0.7 
mm. times. 0.8 mm. times. 1.5 mm along the RL . times .AP . times . SI directions, 
respectively. It took 5 seconds to acquire a 2D TOF image of each axial section. 

Detailed Description Text (45) : 

The k-space sampling pattern for the Gd-enhanced mrMRA was planned in an on-line 
fashion with the help of a scout image of the patient's anatomy. The scout image was 
the axial MIP ("collapsed" image ) of a 2D TOF MRA study of 10 thick axial sections 
covering the carotid bifurcation. It took 50 seconds to acquire these images with 
high resolution (256 . times . 192 ) in the transverse plane. However, an adequate scout 
could easily have been obtained in much less time by increasing slice thickness and 
decreasing resolution. The scout image was used to create a cell array which could 
be translated and sheared so that the arteries were contained in just two cells, 
separated by a single "empty" cell. FIG. 6a shows the resulting cells, which had 
dimensions equal to 24 mm. times. 3 0 mm and belonged to an 8. times. 4 array that 
covered a 190 mm. times. 120 mm region of the FOV. FIG. 6a is an axial* "scout" image 
of the neck of a patient volunteer with a superposed array of cells, just two of 
which contained the carotid and vertebral arteries. FIG. 6b is the sparse pattern of 
862 k-space locations that are optimal for reconstructing an image with edges 
confined to the two vessel -containing cells in panel a. FIG. 6c is a typical 
Gd-enhanced axial image reconstructed from data at the sparse k-space locations in 
panel b. This is the average of two mrMR images , produced from the two 3D mrMRA 
scans that coincided with the first arterial pass of the contrast agent. FIG. 6d are 
MIP views of the carotid bifurcations, created from the complete set of mrMRA images 
like the one in panel c. 

Detailed Description Text (46) : 

FIG. 6e is the corresponding MIPs, reconstructed from unenhanced 2D TOF MRA images . 
Note the loss of TOF flow signal in both carotid bulbs (arrows) which is not 
observed in the Gd-enhanced MIPs. 

Detailed Description Text (47) : 

The relative positions and dimensions of the cells of FIG. 6a were used to prescribe 
t ^ ie k - space sampling pattern to be used in the subsequent Gd-enhanced mrMRA pulse 
sequence. Because the edge component of the image was expected to vanish (or nearly 
vanish) except in the two artery-containing cells, Eqs . (9) and (10) could be used 
to reconstruct it from data on two coarse grids. As in the phantom study, the 
displacement vector between the centers of the two unsheared supporting cells was 
d=(2 . DELTA. x, 0). Therefore, the mrMR sampling theorem dictated that the two coarse 
grids be superposed with the optimal offsets (Eq. (19)): k.sub.l =(1/4 
. DELTA. k. sub. x, 0) and k.sub.2 =0, where . DELTA . k . sub . x =2 .pi . / . DELTA. x. The 
dimensions of these grids were chosen to be 20. times. 16, so that the edge image 
would have spatial resolution equal to 1.2 mm.times.1.9 mm. The image * s background 
component was reconstructed from the data on a dense 24. times. 10 grid at the center 
°f k-space, after multiplying those data by a Hamming filter covering that region. 
FIG. 6b shows the composite mrMR sampling pattern, which was comprised of the two 
large coarse grids, together with the small dense grid. This sampling pattern 
contained just 862 of the 10,240 WKS measurements required to achieve uniform 1.2 
mm.times.1.9 mm resolution throughout the entire FOV. 

Detailed Description Text (48) : 

The patient was infused with 0.3 mM/kg of gadodiamide contrast agent (Omniscan, 
Nycomed, Wayne, Pa.) over a 25 second time period. The mrMRA scan was performed 
repetitively before, during, and after the injection, using a 3D spoiled gradient 
echo pulse sequence (17) with TR=8 ms, TE=1.7 ms, and flip angle=45 . degree . . A 120 
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mm thick coronal slab was selectively excited. It was frequency -encoded along the SI 

direction in order to produce 1.2 mm thick axial sections. The signals were 

phase -encoded along the x (RL) and y (AP) axes according to the composite sampling 

pattern shown in FIG. 6b. It took just seven seconds to acquire the 862 

phase -encoded signals in this pattern. The mrMRA pulse sequence was run seven times 

in quick succession over a 49 second time period that began seven seconds before the 

infusion and ended seventeen seconds after the completion of the infusion. In this 

way, acquisition of a complete mrMR data set was ensured during the first arterial 

pass of contrast agent without the need for "hit or miss" bolus timing techniques. 

Detailed Description Text (49) : 

The signals from non-vascular tissues were suppressed by subtracting the 
pre-infusion mrMRA data from each subsequent data set . In order to identify the 
arterial phase, low resolution images were reconstructed from each subtracted mrMRA 
data set by 2D FT of the signals in the center of k- space . As shown in FIG. 7, these 
images clearly demonstrate early arterial enhancement, followed by venous 
enhancement . 

Detailed Description Text (50) : 

FIG. 7 shows low resolution axial images of the neck of the patient volunteer in 
FIG. 6, at six time intervals after the Gd infusion. These images were produced by 
2D FT of the signals acquired in the center of k-space during each of the last six 
Gd-enhanced mrMRA scans, after the corresponding signals from the first scan were 
subtracted in order to suppress the intensity of stationary tissues. 

Detailed Description Text (51) : 

In order to maximize contrast-to-noise, the subtracted data from the fourth and 
fifth scans was averaged, both of which showed arterial enhancement without 
significant venous signal. The subtracted data in the densely sampled square at the 
center of k-space were - multiplied by the Hanmming filter that covered that area, and 
the background component of the mrMR image was reconstructed by FFT of this filtered 
data . The edge component of the mrMR image was reconstructed from the subtracted 
data on the two coarse k-space grids, after the data had been multiplied by the 
complement of the Hamming filter (see Eq. (15)) and Fourier transformed. The 
complete mrMR image was calculated by adding the resulting edge component to the 
background component (FIG. 6c) . Finally, MIPs of the mrMRA reconstructions were 
created and compared to MIPs of the unenhanced 2D TOF images acquired prior to 
infusion (FIGS. 6d and e) . 

Detailed Description Text (54) : 

FIG. le shows the mrMR reconstruction of the simulated image in FIG. la. The two 
images are identical to within round-off error, despite the fact that the former 
image was reconstructed from 1/12 of the data used to reconstruct the latter one. 
FIG. 2g shows the mrMR reconstruction of the simulated image in FIG. 2a. Notice the 
accuracy of the result, even though it was reconstructed from sparse data (FIGS. 2e 
and 2f) which comprised only 1/9 of the complete FT data set . FIG. 2h is the error 
image, which is the magnitude of the difference between the mrMR image and the exact 
image . This image is windowed to show the discrepancies. These errors were caused by 
some low spatial frequencies in the edge component which were outside the three 
supporting cells and led to aliasing when this edge image was reconstructed by mrMR. 
Notice that this error is comprised solely of low frequency image components; i.e. 
the image ! s edges are reconstructed faithfully. 

Detailed Description Text (56) : 

FIG. 5a shows the exact image of a typical slice at the level of the tubing 
bifurcations, reconstructed from the complete set of 16,384 data required to image 
the entire FOV with an isotropic 1.25 mm resolution. FIG. 5c shows the mrMR i?}^.9^. 
that was created from sparse data at the k'.sgajse locations in FIG. 5b, which 
comprised less than 1/8 of the complete^ data set. FIG. 5d is the error in?s.9.€J/„ 
equal to the magnitude of the difference between the mrMR Aj:A a i3.?. anc ^ the exac £ ^{??.9.?..- 
This image is windowed to increase the conspicuity of the faint 

The error image also shows truncation s^.ti^acts. at the surface of the bath, which 
was outside the two edcje -containing cells. FIGS. 5e and 5f show maximum intensity 
projections (MIPs) created from the exact and mrMR reconstructions, respectively. 
Notice that the mrMR MIP accurately depicts the tiny " .rite^_-of f 11 discontinuities of 
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the lumen where the caliber of the tubing changes abruptly. These features resemble 
minimal atherosclerotic plaques. Only faint artifacts are present, as demonstrated 
by the MIP error image (FIG. 5g) , which is windowed to accentuate these 
discrepancies . 

Detailed Description Text (58) : 

FIG. 6c shows a typical axial mrMR image that was reconstructed from the data at the 
862 k- space locations in FIG. 6b. FIG. 6d shows MIP images derived from mrMR images 
of 60 such axial sections. These should be compared to the MIPs of the unenhanced 2D 
TOF images in FIG. 6e. Notice that the 2D TOF MIPs are degraded by signal loss in 
the carotid bulbs, which are the most common locations of atherosclerotic plaques. 
In contrast, the mrMRA MIPs depict the full extent of the carotid bulbs. 

Detailed Description Text (60) : 

If the FT of a signal is limited to a single band, the WKS sampling theorem shows 
how to collect the minimum number of samples from which the entire signal can be 
reconstructed without noise amplification. Conventional FT MR I is based on this 
theorem. A more general sampling theorem is described below which shows how to 
sample the k- space of an image that vanishes except in multiple regions (i.e. it is 
"multi -banded" in image space ) . These data can be used to reconstruct the image 
without noise amplification; i.e. the image has the same noise level as a comparable 
conventional FT image with the same spatial resolution over the same total area of 
support. Furthermore, it is demonstrated below that this method can be modified to 
handle a larger class of images : those that have slowly varying signal intensity 
across the entire FOV and high-contrast edges confined to small regions. The 
"background" component of such an image can be reconstructed from densely sampled 
data in the center of k- space, while its "edge " component can be reconstructed from 
the sparse k- space samples dictated by the new sampling theorem. 

Detailed Description Text (61) : 

If the edge -containing regions cover only a small portion of the FOV, the image can 
be reconstructed without noise amplification from a sparse subset of the complete 
(WKS) data set . This makes it possible to use prior knowledge of the object to 
increase the efficiency of the sampling process . This can be advantageous in several 
ways. In many cases, mrMR images can be acquired in less time (or with greater 
spatial resolution) because it takes fewer TR intervals to sample the image at the 
sparse mrMR k- space locations than at the more numerous WKS points. Furthermore, if 
k- space trajectories are used to collect multiple k - space samples per TR interval, 
lower bandwidth and/or shorter read-out time (trajectory traversal time) can be used 
because it is only necessary to visit the sparse k - space locations needed by mrMRI . 
Reducing the bandwidth leads to increased signal-to-noise, while shortening the 
read-out time tends to reduce artifacts due to of f -resonance effects. 

Detailed Description Text (62) : 

It should be noted that the mrMR method is computationally cheap, because images are 
formed from linear combinations of small FFTs of the data. F ( or imaging situations in 
which the supporting cells occupy only a small fraction of the entire FOV , the 
reconstruction process is computationally cheaper than that of conventional FT MRI , 
because the latter involves a single large FFF. In the experiments discussed here, a 
single R10000 processor 104 (Silicon Graphics, Inc., Mountain View, Calif.) computed 
each mrMR image in a small fraction of a second. In fact, in the current 
implementation, the speed of mrMR image reconstruction was limited by extraneous 
factors, such as the time required for the Unix script to call the reconstruction 
routine and the time required to print a message on the user's screen. 

Detailed Description Text (63) : 

The mrMR method was applied to simulated data and to experimental data from a 
vascular phantom. When the resulting .±-™&3®&„ were compared to "exact" iroag_es,_ it was 
apparent that there was no noise amplification and that little .?AA a A?.il L 3.. or 
truncation &£tif.act was produced by the presence of a slowly undulating background 
intensity. It was also shown how mrMRI could be used to perform first -pass 
Gd-enhanced 3D carotid MRA in just 7 seconds. The exam was performed in a fashion 
analogous to conventional angiography: namely, a series of 3D mrMRA scan_s was 
performed before, during, and after the injection of contrast agent. The resulting 
data were used to retrospectively identify the first arterial pass of contrast. In 
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other words, it was not necessary to "catch" the arterial phase with bolus timing 
procedures, which can be operator-dependent and patient -dependent (13-16) . Bolus 
timing can be particularly ticklish for Gd-enhanced carotid imaging because the 
internal jugular veins enhance within 10 seconds of the first arterial pass of 
contrast agent and because they are immediately adjacent to the carotid arteries. 

Detailed Description Text (64) : 

The mrMRI method may also be useful in other imaging situations in which the high 
spatial frequencies are known to be confined to multiple regions of the FOV . For 
example, this type of prior knowledge may be available in functional MRI, 
interventional MRI, and MR spectroscopic imaging, as well as in other applications 
of MRA. Three-dimensional MR spectroscopic imaging may benefit from extensions of 
mrMRI to higher dimensions. Moreover, the new sampling theorems might increase the 
sampling efficiency of other Fourier imaging modalities, such as radio astronomy, 
radar, and crystallography. It may also be applied in telecommunications for the 
purpose of digitizing and reconstructing signals with frequency components in 
multiple frequency bands. 

Detailed Description Text (65) : 

There are two principle sources of image degradation in mrMRI reconstructions. First 
of all, the SNR of mrMR images is limited by the total number of measured signals, 
just as in conventional FT reconstruction. The SNR of our carotid mrMRA images was 
suboptimal because they were reconstructed from just 862 k-space samples and because 
no attempt was made to optimize our pulse sequence and infusion protocol. Secondly, 
any high spatial frequency features (i.e. edges ) which fall outside of our specified 
cells (thereby violating our assumptions) will alias into these cells. These 
artifacts can be minimized by using techniques that increase the contrast of edges 
within the specified cells and simultaneously suppress the contrast of edges outside 
of these cells. Image quality may also be improved by: 1) optimization of the 
contrast-to-noise of the pulse sequence (e.g. use of fat suppression, fractional 
echoes, and more suitable values of the TR, bandwidth, and flip angle parameters) ; 
2) utilization of a power injector in order to increase the concentration of 
gadodiamide in the blood; 3) use of better spatial filters for separating the image 
into background and edge components; 4) creation of the mrMRA images from 
temporally-filtered signals collected during all seven mrMR scans (not just during 
one or two arterial acquisitions as in IIIC) . By making use of data from all of the 
mrMRA acquisitions, significantly higher SNR may be achieved (e.g., 2-3 times 
higher) . Similar temporal filtering methods have been developed for conventional 
X-ray angiography. 

Detailed Description Text (67) : 

It is instructive to compare mrMRI to Locally Focused MR I (LF MRI) , also developed 
by the inventors. In most situations, the mrMR method is preferable because it 
avoids the noise amplification of LF MRI and because it is computationally much 
cheaper than LF MRI. In order to perform LF MRI, the first step is to generate a set 
of non-Fourier basis functions that span the space of all images having the pattern 
of spatial variation suggested by prior knowledge. The Fourier components of the 
image are linearly related to the image 1 s projections onto these LF basis functions. 
Therefore, these linear equations can be inverted to calculate the basis function 
coefficients in terms of an equal number of measured Fourier components. If our 
prior knowledge is strong enough to be expressed by a small number of basis 
functions, the image can be reconstructed from a small number of k •- space samples. 
Notice that the LF MRI method is quite general, allowing the user to exploit prior 
knowledge of the image 1 s general pattern of spatial variation in multiple regions of 
irregular shape. The main problem is that LF images are degraded by noise 
amplification, caused by ill -conditioning of the matrix inversion required to find 
the basis function coefficients. The elements of this matrix depend on the !$" 8 .E?. C .? 
sampling pattern, as well as on the form of the basis functions dictated by our 
prior knowledge. In general, it has not been possible to find minimal 
sampling patterns that make this matrix well-conditioned (i.e. that make its 
singular values approximately equal) . In order to make the problem better 
conditioned, it has been necessary to over -determine its solution by sampling 4-6 
times more k-space points than the number of basis functions. Of course, this is 
undesirable because it increases scan time by a factor of 4-6 above the minimal LF 
MRI requirements. The second disadvantage of the LF technique is its computational 
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expense. Iruage reconstruction requires the inversion of a matrix that has dimensions 
equal to the number of basis functions, which is itself proportional to the total 
area of the edge -containing regions. For example, a 400 . times . 400 matrix was 
inverted in order to reconstruct Gd-enhanced MRA images of the carotid arteries with 
isotropic 1.25 mm resolution. It took approximately 3 0 seconds for a single R10000 
processor (Silicon Graphics, Inc., Mountain View, Calif.) to perform this 
computation with a QR factorization algorithm. In many applications, this 
computation may have to be performed only once for a given patient; i.e. it need not 
be repeated on a slice-by-slice basis if the k-space sampling and prior knowledge 
are the same for all slices. However, the expense of this computation increases 
rapidly with the matrix dimension. Therefore, it could become a limiting factor if 
the prior knowledge is not strong enough to effectively limit the number of basis 
functions . 

Detailed Description Text (68) : 

Notice that mrMRI utilizes prior knowledge of the approximate locations of edges , 
but no assumptions are made about the configurations of those edges . Therefore, 
mrMRI is less model -dependent than methods "which assume that the edge configurations 
of the unknown image are similar to those of prior images of the subject. For 
example, in one prior publication entitled "Reduced encoding Imaging by Generalized 
series Reconstruction" (RIGR) , the image is expressed in terms of a small number of 
basis functions that are products of slowly varying Fourier functions and a high 
resolution baseline image of the subject's anatomy. Therefore, all edges in a RIGR 
image must have the same configuration as the edges in the anatomical reference 
image , and problems can be expected if the unknown image contains new sharp edges . 
In "Feature-Recognizing MRI", a reduced basis set is identified as the principal 
components of "training" images of other subjects, who are thought to resemble the 
unknown subject. Images of the unknown subject are found by computing the basis 
function' coefficients from a reduced set of k-space samples. Of course, these images 
will only be accurate to the extent to which the unknown subject resembles the 
collection of training subjects. Furthermore, both "feature-recognizing" and RIGR 
images will suffer from noise amplification because the computation of basis 
function coefficients is usually not well -conditioned . In "Singular Value 
Decomposition MRI" (SVD MRI), a reduced basis set is derived from the singular value 
decomposition (SVD) of a single prior image of the subject. This is mathematically 
equivalent to applying the "feature -recognizing" procedure to the collection of 
columns in that image . Other investigators have found that this method is 
model -dependent in the following sense: the SVD basis functions may not accurately 
represent new edges which are present in the image of interest but were not present 
in the prior image . Furthermore, SVD MRI also suffers from inherently low SNR. 
Zientara et al implemented SVD MRI by using direct RF excitation to encode signals 
with the profiles of the basis functions. In this way, the object's projection onto 
each basis function can be measured directly, and it is possible to avoid the 
ill-conditioned computation of the basis function coefficients from Fourier encoded 
signals. Unfortunately, direct excitation requires the use of variable flip angles 
across the excitation profile, and this may produce 3-4 times less signal than that 
obtained with a uniform flip angle excitation followed by Fourier encoding . 

Detailed Description Text (70) : 

In another illustrated embodiment, the method of reconstructing images is extended 
to time dependent signals. As used herein, a time dependent signal may be any 
information-containing signal transmitted in the radio frequency, microwave, 
infrared, or visible range. 

Detailed Description Text (71) : 

Tne process of reconstructing time dependent signals is very similar to 
reconstruction of MR and other image s . For example, MR images are reconstructed from 
samples of the image's Fourier representation at certain locations in the Fourier 

0<- space ) . The frequency representation of a time dependent signal can be 
considered to be a one -dimensional image . The Fourier representation of that ID 
image is the time dependent signal itself. Just as in MRI, such a ID image can be 
reconstructed from samples of its Fourier representation at certain locations in the 
Fourier space; in other words, the signal's frequency representation can be 
reconstructed from samples of the time dependent signal at certain times. Then, the 
resulting frequency representation of the signal can be Fourier transformed to 
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reconstruct the time -dependent signal at any time. 
Detailed Description Text (72) : 

The method of sampling and reconstructing time dependent signals may be used to 
transmit, receive, and reconstruct signals that have frequency components in 
multiple, possibly non-adjacent frequency bands. For example, it may be desirable to 
transmit a certain amount of information to a receiver over a certain period of 
time. However, there may not be an available band of frequency that is broad enough 
to permit the information to be transmitted in the desired time interval. If a 
sufficient number of narrow bands of frequency are available, the described 
invention makes it possible to use those multiple narrow frequency bands to transmit 
the information in the desired time interval. This process is described in the 
following example. 

Detailed Description Text (76) : 

The receiving unit would use the known locations and breadths of the narrow 
frequency bands to determine the times at which the filtered signal should be 
sampled digitally. In general, the temporal sampling pattern would be comprised of a 
superposition of coarse sampling patterns. Each such coarse pattern would be 
suitable for Fourier transform reconstruction of any signal with frequency 
components in just one of the narrow bands. Also, each such coarse pattern would 
have a different temporal off -set, such off -sets being chosen in order to make the 
eigenvalues of the reconstruction matrix comparable to one another. 

Detailed Description Text (78) : 

The receiving unit would Fourier transform subsets of the stored sampled data and 
use the reconstruction matrix to linearly combine those Fourier transforms in order 
to reconstruct the signal's frequency components in each narrow frequency band. 

Detailed Description Text (80) : 

A specific embodiment of a novel method and apparatus for reconstructing images 
according to the present invention has been described for the purpose of 
illustrating the manner in which the invention is made and used. It should be 
understood that the implementation of other variations and modifications of the 
invention and its various aspects will be apparent to one skilled in the art, and 
that the invention is not limited by the specific embodiments described. Therefore, 
it is contemplated to cover the present invention any and all modifications, 
variations, or equivalents that fall within the true spirit and scope of the basic 
underlying principles disclosed and claimed herein. 

Other Reference Publication (1) : 

Beaty, M.G., et al . , Derivative Sampling for Multiband Signals, Numer. Funct . Anal, 
and Qptimiz . , 20(9&10), 875-898 (1989). 

Other Reference Publication (3) : 

Beatty, M.G., Multichannel sampling for multiband signals, 36 Signal Processing 
133-138 (1994) . 

Other Reference Publication (13) : 

Walter, G.G., Recent Extension of the Sampling Theorem, Signal Processing, Part I: 
Signal Processing Theory (Auslander et al . , eds.) 229-238 (Springer-Verlag, 1990). 



1. A method of acquiring and reconstructing an image, such method comprising the 
steps of: 

obtaining prior knowledge of the .image, possibly by coarse sampling of the -"i-"iage_; 
using the obtained prior knowledge of the image to identify relative locations of 
structures having relatively high contrast > 

prescribing a' set oil k-sgace locations based upon the relative locations of the 
structures in order to achieve comparable eigenvalues of a reconstruction matrix; 

sampling the k- space at the prescribed k- space locations to obtain k- space sample 
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data; 



decomposing the k- space sample data into background data and edge data; 



Fourier transforming the background data to reconstruct a background 
component ; 



image 



Fourier transforming subsets of the edge data and using the reconstruction matrix to 
form a linear combination of these Fourier transformations in order to reconstruct 
an edge image component; and 

combining the background image component and the edge image component to generate a 
final image . 

2 . The method of acquiring and reconstructing an image as in claim 1 further 
comprising selecting a grid-like pattern of cells which substantially segregates 
high contrast areas, having the high contrast edges, from other areas. 

3 . The method of acquiring and reconstructing an image as in claim 2 further 
comprising prescribing a plurality of k-space sampling patterns for the high 
contrast areas and for the other areas, one of such k-space sampling patterns 
comprising a dense collection of k-space locations near the center of k-space, and 
each of the other such k-space sampling patterns comprising a coarse grid of k-space 
locations covering a relatively large region of k-space, each such coarse grid being 
suitable for Fourier transform reconstruction of any image having features in just 
one cell and each such coarse grid having a pre-determined off -set in k-space . 

4 . The method of acquiring and reconstructing an image as in claim 3 further 
comprising determining the reconstruction matrix based upon the locations of the 
segregated high contrast areas and the off -sets of the prescribed coarse k-space 
sampling patterns. 

5. The method of acquiring and reconstructing an image as in claim 4 further 
comprising computing the eigenvalues of the determined reconstruction matrix. 

6. The method of acquiring and reconstructing an image as in claim 5 further 
comprising adjusting the off -sets of the prescribed coarse k- space sampling patterns 
so that the computed eigenvalues of the determined reconstruction matrix are 
approximately equal . 

7. The method of acquiring and reconstructing an image as in claim 6 further 
comprising sampling the data at all of the locations in all of the prescribed dense 
and coarse k- space sampling patterns. 

8. The method of acquiring and reconstructing an image as in claim 7 further 
comprising multiplying data at the locations of the dense k- space pattern by values 
corresponding to a low-pass filter to provide low-pass filtered data. 

9. The method of acquiring and reconstructing an image as in claim 7 further 
comprising multiplying data at the locations of the coarse k-space sampling patterns 
by values corresponding to a complementary high-pass filter to provide high-pass 
filtered data. 

10. The method of acquiring and reconstructing an image as in claim 8 further 
comprising Fourier transforming the low-pass filtered data to reconstruct the 
background image component . 

11. The method of acquiring and reconstructing an ; }„ m .^.9?. as in cla i m 9 further 
comprising Fourier transforming the high-pass filtered data of each coarse k-space 
sampling pattern and using the reconstruction matrix to linearly combine the Fourier 
transforms to reconstruct the edge component of the image- 

12. The method of acquiring and reconstructing an image as in claim 11 further 
comprising adding the background image component to the edoe „„imaoe component to 
reconstruct the i mage . 
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13. A method of digitizing and reconstructing a time -dependent signal, such method 
comprising the sbeps of: 

obtaining prior knowledge of frequencies at which the signal has Fourier components 
which exceed a threshold, possibly by filtering the signal to remove or reduce all 
other undesired frequency components; 

using the obtained prior knowledge to prescribe a set of times in order to achieve 
comparable eigenvalues of a reconstruction matrix; 

sampling the signal at the prescribed times in order to obtain data; 

Fourier transforming subsets of the sampled data and using the reconstruction matrix 
to form a linear combination of these Fourier transformations in order to 
reconstruct the signal's Fourier components in each frequency band; and 

inverse Fourier transforming the reconstructed Fourier components in each frequency 
band in order to reconstruct the time -dependent signal having frequencies in each 
frequency band. 

15. The method of digitizing and reconstructing a time -dependent signal as in claim 

14 further comprising prescribing a plurality of temporal sampling patterns for the 
spectral bands having frequency activity, each such temporal pattern consisting of a 
coarse pattern of times which pattern is suitable for Fourier transform 
reconstruction of the frequency components of any signal having just one frequency 
band and which pattern has a predetermined off -set in time. 

16. The method of digitizing and reconstructing a time-dependent signal as in claim 

15 further comprising determining the reconstruction matrix depending on the 
frequency locations of the segregated spectral bands and the time off -sets of the 
prescribed coarse temporal sampling patterns. 

17. The method of digitizing and reconstructing a time-dependent signal as in claim 

16 further comprising computing the eigenvalues of the determined reconstruction 
matrix. 

18. The method of digitizing and reconstructing a time -dependent signal as in claim 

17 further comprising adjusting the prescribed time off -sets so that the computed 
eigenvalues of the determined reconstruction matrix are approximately equal . 

22 • Apparatus for acquiring and reconstructing an image, such apparatus comprising: 

means for obtaining prior knowledge of the image, possibly by coarse sampling of the 
image ; using the obtained prior knowledge of the image to identify relative 
locations of structures having relatively high contrast edges ; 

means for prescribing a set of k- space locations based upon the relative locations 
of the structures in order to achieve comparable eigenvalues of a reconstruction 
matrix; 

means for sampling the k- space at the prescribed k- space locations to obtain k- space 
sample data; 

means for decomposing the k- space sample data into background data and edge data; 

means for Fourier transforming the background data to reconstruct a background Ml a H.5L 
component ; 

means for Fourier transforming subsets ^ data and using the reconstruction 

matrix to form a linear combination of these Fourier transformations in order to 
reconstruct an edge ifnage. component; and 

means for combining the background component and the .?49.5..All!?E5.. component to 

generate a final ,i^S. e ..- 
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23 . The apparatus for acquiring and reconstructing an image as in claim 22 further 
comprising means for selecting a grid-like pattern of cells which substantially 
segregates high contrast areas, having the high contrast edges, from other areas. 

24 . The apparatus for acquiring and reconstructing an image as in claim 23 further 
comprising means for prescribing a plurality of k- space sampling patterns for the 
high contrast areas and for the other areas, one of such k- space sampling patterns 
comprising a dense collection of k- space locations near the center of k- space, and 
each of the other such k- space sampling patterns comprising a coarse grid of k-space 
locations covering a relatively large region of k-space, each such coarse grid being 
suitable for Fourier transform reconstruction of any image having features in just 
one cell and each such coarse grid having a pre-determined off -set in k-space . 

25. The apparatus for acquiring and reconstructing an image as in claim 24 further 
comprising means for determining the reconstruction matrix based upon the locations 
of the segregated high contrast areas and the off -sets of the prescribed coarse 
k-space sampling patterns. 

26. The apparatus for acquiring and reconstructing an image as in claim 25 further 
comprising means for computing the eigenvalues of the determined reconstruction 
matrix. 

27. The apparatus for acquiring and reconstructing an image as in claim 26 further 
comprising means for adjusting the off -sets of the prescribed coarse k-space 
sampling patterns so that the computed eigenvalues of the determined reconstruction 
matrix are approximately equal . 

28. The apparatus for acquiring and reconstructing an image as in claim 27 further 
comprising means for sampling the data at all of the locations in all of the 
prescribed dense and coarse k-space sampling patterns. 

29. The apparatus for acquiring and reconstructing an image as in claim 28 further 
comprising means for multiplying data at the locations of the dense k-space pattern 
by values corresponding to a low-pass filter to provide low-pass filtered data. 

30. The apparatus for acquiring and reconstructing an image as in claim 29 further 
comprising means for multiplying data at the locations of the coarse k-space 
sampling patterns by values corresponding to a complementary high-pass filter to 
provide high-pass filtered data. 

31. The apparatus for acquiring and reconstructing an image as in claim 29 further 
comprising means for Fourier transforming the low-pass filtered data to reconstruct 
the background image component . 

32. The apparatus for acquiring and reconstructing an image as in claim 30 further 
comprising means for Fourier transforming the high-pass filtered data of each coarse 
k- space sampling pattern and using the reconstruction matrix to linearly combine the 
Fourier transforms to reconstruct the edge component of the image . 

33. The apparatus for acquiring and reconstructing an image as in claim 32 further 
comprising means for adding the background image component to the edge image 
component to reconstruct the image . 

34 * Apparatus for digitizing and reconstructing a time-dependent signal, such 
apparatus comprising: 

means for obtaining prior knowledge of the frequencies at which the signal has 
significantly large Fourier components, possibly by filtering the signal to remove 
or reduce all other frequency components; 

means for using the obtained prior knowledge to prescribe a set of times in order to 
achieve comparable eigenvalues of a reconstruction matrix; 

means for sampling the signal at the prescribed times in order to obtain data; 
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means for Fourier transforming subsets of the sampled data and using the 
reconstruction matrix to form a linear combination of these Fourier transformations 
in order to reconstruct the signal's Fourier components in each frequency band; and 

means for inverse Fourier transforming the reconstructed Fourier components in each 
frequency band in order to reconstruct the time-dependent signals with frequencies 
in each frequency band. 

35. The apparatus for digitizing and reconstructing a time-dependent signal as in 
claim 34 further comprising means for selecting a spectral pattern which 
substantially segregates spectral bands having frequency activity above a 
predetermined threshold, from other areas. 

36. The apparatus for digitizing and reconstructing a time -dependent signal as in 
claim 35 further comprising means for prescribing a plurality of temporal sampling 
patterns for the spectral bands having frequency activity, each such temporal 
pattern consisting of a coarse pattern of times which pattern is suitable for 
Fourier transform reconstruction of the frequency components of any signal having 
just one frequency band and which pattern has a predetermined off -set in time. 

37. The apparatus for digitizing and reconstructing a time-dependent signal as in 
claim 3 6 further comprising means for determining the reconstruction matrix 
depending on the frequency locations of the segregated spectral bands and the time 
of f -sets of the prescribed coarse temporal sampling patterns. 

38. The apparatus for digitizing and reconstructing a time -dependent signal as in 
claim 37 further comprising means for computing the eigenvalues of the determined 
reconstruction matrix. 

39. The apparatus for digitizing and reconstructing a time-dependent signal as in 
claim 3 8 further comprising means for adjusting the prescribed time off -sets so that 
the computed eigenvalues of the determined reconstruction matrix are approximately 
equal . 

40. The apparatus for digitizing and reconstructing a time-dependent signal as in 
claim 3 9 further comprising means for using an analogue or digital filter to 
suppress signals with frequency components outside the segregated spectral bands and 
for sampling the filtered signal at all locations in all of prescribed coarse 
temporal sampling patterns and for storing the sampled data. 

41. The apparatus for digitizing and reconstructing a time-dependent signal as in 
claim 4 0 further comprising means for Fourier transforming the sampled data in each 
coarse temporal sampling pattern and using the reconstruction matrix to linearly 
combine the Fourier transforms to reconstruct the frequency components in each 
frequency band. 

42. The apparatus for digitizing and reconstructing a time-dependent signal as in 
claim 41 further comprising means for inverse Fourier transforming the reconstructed 
frequency components in each frequency band in order to reconstruct the 
time-dependent signals having frequencies in each frequency band. 

43 * Apparatus for digitizing and reconstructing a time-dependent signal, such 
apparatus comprising: 

a scanning receiver which obtains prior knowledge of the frequencies at which the 
signal has significantly large Fourier components, possibly by filtering the signal 
to remove or reduce all other frequency components; 

a matrix g£Ocessor_ which uses the obtained prior knowledge to prescribe a .set of 
times in order to achieve comparable eigenvalues of a reconstruction matrix; 

an analogue or digital filtering P£°c.essor__ which suppresses signals with frequency 
components outside the spectral bands of interest; 



19 of 20 



10/27/2003 9:01AM 



Record Display Form 



http://westbrs: 8002/bin/gate . exe?f=do^^ . . c_2=&p_doc_3 =&p_doc_4=&p_doc__5 =&p_doc_6= 




a sampling processor' which samples the filtered signal at the prescribed times in 
order to obtain data and which stores the sampled data; 

a Fourier processor which Fourier transforms subsets of the sampled data and uses 
the reconstruction matrix to form a linear combination of these Fourier 
transformations in order to reconstruct the signal's Fourier components in each 
frequency band; and 

an inverse Fourier processor which inverse Fourier transforms the reconstructed 
Fourier components in each frequency band in order to reconstruct the time -dependent 
signals in each frequency band. 

44. The apparatus for digitizing and reconstructing a time -dependent signal as in 
claim 43 further comprising a wide-band receiver. 

45. The apparatus for digitizing and reconstructing a time-dependent signal as in 
claim 44 further comprising an analog or digital filter which suppresses radio 
frequency energy outside desired spectral ranges. 

46. The apparatus for digitizing and reconstructing a time-dependent signal as in 
claim 45 further comprising a sampling device which samples the filtered signal at 
the prescribed times, and stores the resulting data. 

47. The apparatus for digitizing and reconstructing a time-dependent signal as in 
claim 46 further comprising a reconstruction processor which reconstructs the 
time-dependent signal, which has frequency components in each frequency band, from 
the sampled data. 
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Jul 27, 1999 



DOCUMENT- IDENTIFIER: US 5928148 A 

TITLE: Method for performing magnetic resonance angiography over a large field of 
view using table stepping 

Abstract Text (1) : 

MRA data is acquired from a large region of interest by translating the patient to 
successive stations at which successive portions of the MRA data set are acquired. 
Patient movement is chosen to track a bolus of contrast agent as it passes through 
the region of interest to achieve maximum image contrast. In one embodiment a 
stationary local coil is supported adjacent the patient to acquire the MRA data and 
in another embodiment a multi -segment local coil moves with the patient and its 
segments are sequentially switched into operation. 

Parent Case Text (2) : 

This application is based on provisional application Ser. No. 60/048,286, filed on 
Jun. 2, 1997 and entitled " METHOD FOR PERFORMING MAGNETIC RESONANCE ANGIOGRAPHY OVER 
A LARGE FIELD OF VIEW USING TABLE STEPPING " . 

Brief Summary Text (3 ) : 

Diagnostic studies of the human vasculature have many medical applications. X-ray 
imaging methods such as digital subtraction angiography ("DSA") have found wide use 
in the visualization of the cardiovascular system, including the heart and 
associated blood vessels. One of the advantages of these x-ray techniques is that 
image data can be acquired at a high rate (i.e. high temporal resolution) so that a 
sequence of images may be acquired during injection of the contrast agent. Such 
"dynamic studies" enable one to select the image in which the bolus of contrast 
agent is flowing through the vasculature of interest. Images showing the circulation 
of blood in the arteries and veins of the kidneys, the neck and head, the 
extremities and other organs have immense diagnostic utility. Unfortunately, 
however, these x-ray methods subject the patient to potentially harmful ionizing 
radiation and often require the use of an invasive catheter to inject a contrast 
agent into the vasculature to be imaged / There is also the issue of increased 
nephro-toxicity and allergic reactions to iodinated contrast agents used in 
conventional x-ray angiography. 

Brief Summary Text (4) : 

Magnetic resonance angiography (MRA) uses the nuclear magnetic resonance (NMR) 
phenomenon to produce images of the human vasculature . When a substance such as 
human tissue is subjected to a uniform magnetic field (polarizing field B.sub.O), 
the individual magnetic moments of the spins in the tissue attempt to align with 
this polarizing field, but precess about it in random order at their characteristic 
Larmor frequency. If the substance, or tissue, is subjected to a magnetic field 
(excitation field B.sub.l) which is in the x-y plane and which is near the Larmor 
frequency, the net aligned moment, M.sub.z, may be rotated, or "tipped", into the 
x-y plane to produce a net transverse magnetic moment M.sub.t. A signal is emitted 
by the excited spins, and after the excitation signal B.sub.l is terminated, this 
signal may be received and processed to form an image . 

Brief Summary Text (5) : 

When utilizing these signals to produce images, magnetic field gradients (G.sub.x 
G.sub.y and G.sub.z) are employed. Typically, the region to be imaged is scanned by 
a sequence of measurement cycles in which these gradients vary according to the 
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particular localization method being used. The resulting set of received NMR signals 
are digitized and processed to reconstruct the image using one of many well known 
reconstruction techniques. 

Brief Summary Text (6) : 

MR angiography (MRA) has been an active area of research. Two basic techniques have 
been proposed and evaluated. The first class, time-of -flight (TOF) techniques, 
consists of methods which use the motion of the blood relative to the surrounding 
tissue. The most common approach is to exploit the differences in magnetization 
saturation that exist between flowing blood and stationary tissue. Flowing blood, 
which is moving through the excited section, is continually refreshed by spins 
experiencing fewer excitation pulses and is, therefore, less saturated. The result 
is the desired image contrast between the high-signal moving blood and the 
low-signal stationary tissues. 

Brief Summary Text (7) : 

MRA methods have also been developed that encode motion into the phase of the 
acquired signal as disclosed in U.S. Pat. No. Re. 32,701. These form the second 
class of MRA techniques and are known as phase contrast (PC) methods. Currently, 
most PC MRA techniques acquire two images, with each image having a different 
sensitivity to the same velocity component. Angiographic images are then obtained by 
forming either the phase difference or complex difference between the pair of 
velocity -encoded images . 

Brief Summary Text (8) : 

To enhance the diagnostic capability of MRA a contrast agent such as gadolinium can 
be injected into the patient prior to the MRA scan . Excellent diagnostic images may 
be acquired using contrast -enhanced MRA if the data acquisition is properly timed 
with the bolus passage. 

Brief Summary Text (9) : 

The non- invasiveness of MRA makes it a valuable screening tool for cardiovascular 
diseases. Screening typically requires imaging vessels in a large volume. This is 
particularly true for diseases in the runoff vessels of the lower extremity. The 
field of view ( POV ) in MR imaging is limited by the volume of the B.sub.O field 
homogeneity and the receiver coil size (typically, the F0V<48 cm on current 
commercial MR scanners ) . The anatomic region of interest in the lower extremity, for 
example, is about 100 cm and this requires several POVs , or stations, for a complete 
study. This requires that the patient be repositioned inside the bore of the magnet, 
the patient be relandmarked, scout images be acquired and a preparation scan be 
performed for each FQV . All of' these additional steps take time and, therefore, are 
expensive. When contrast enhanced MRA is performed, the repositioning also 
necessitates additional contrast injections. 

Brief Summary Text (11) : 

The present invention is a method for acquiring NMR data from a large region of 
interest by acquiring NMR data from a series of smaller fields of view which 
collectively span the large region of interest. The patient is automatically 
translated, or stepped, to a new position within the bore of the magnet by moving 
the patient table between the acquisition of each field of view. Scan parameters 
remain constant throughout the procedure and the separate reconstructed images are 
registered and combined to provide a single image of the large region of interest. 

Brief Summary Text (12) : 

Another aspect of the present invention is to acquire the entire region of interest 
with one injection of contrast agent. The acquisition of each field of view is timed 
to correspond with the peak image contrast produced by the bolus of contrast agent . 
The order in which the separate fields of view are acquired and the speed of the 
patient translation between fields of view are selected to track the peak 
intraarterial contrast as the bolus transits through the region of interest. 

Brief Summary Text (13) : 

Yet another aspect of the present invention is the use of local coils to improve the 
signal-to-noise ratio of the acquired MRA image . In one embodiment an array of local 
coils are positioned adjacent the patient and span the large region of interest. 



2 of 10 



10/27/2003 8:59 AM 



Record Display Form http://weslbrs: 8002/biiVgate. exe?f=do^^. . . c_2=&p_doc_3=&p_doc_4-&p_doc_5=&p_doc_6= 



This local coil array moves with the patient as the table is translated from one 
station to the next, and successive segments of the array are connected to the 
receiver as they move into the imaging field of view. In a second embodiment a 
stationary local coil is supported adjacent the patient in the imaging field of 
view. The patient is translated through the imaging field of view and the stationary 
reception field of the local coil. 

Brief Summary Text (14) : 

A general object of the invention is to efficiently acquire MRA data from a large 
region of interest. A single scan is performed in which MRA data is acquired from 
successive parts of the large region of interest. This is accomplished by 
translating the patient between acquisitions and concatenating the acquired data to . 
form a single image of the entire region of interest. 

Brief Summary Text (15) : 

Another object of the invention is to provide optimal contrast by timing the 
acquisition of MRA data with the arrival of contrast agent at each successive FOV in 
the region of interest. This is accomplished by acquiring data from which the 
velocity of the contrast agent bolus as it transits the region of interest is 
estimated. This velocity estimate is used to control the rate at which the patient 
is translated to successive FOV s during the scan and to track the peak in the 
contrast produced by the bolus . 

Brief Summary Text (16) : 

Yet another object of the invention is to improve the SNR of the acquired image . 
This is accomplished by using a local coil which is positioned adjacent the patient. 
In one embodiment the local coil remains stationary as the patient is translated by 
it during the scan . In another embodiment a multi-segment local coil is carried by 
the patient. The local coil segments are successively switched into operation as the 
patient translates through the bore during the scan. 

Drawing Description Text (7) : 

FIG. 6 is a flow chart of the sequence of steps performed by the MRI system of FIG. 
1 to practice the preferred embodiment of the invention; 

Detailed Description Text (2) : 

Referring first to FIG. 1, there is shown the major components of a preferred MRI 
system which incorporates the present invention. The operation of the system is 
controlled from an operator console 100 which includes a keyboard and control panel 
102 and a display 104. The console 100 communicates through a link 116 with a 
separate computer system 107 that enables an operator to control the production and 
display of images on the screen 104. The computer system 107 includes a number of 
modules which communicate with each other through a backplane. These include an 
image processor module 106, a CPU module 108 and a memory module 113, known in the 
art as a frame buffer for storing image data arrays. The computer system 107 is 
linked to a disk storage 111 and a tape drive 112 for storage of image data and 
programs , and it communicates with a separate system control 122 through a high 
speed serial link 115. 

Detailed Description Text (3) : 

The system control 122 includes a set of modules connected together by a backplane. 
These include a CPU module 119 and a pulse generator module 121 which connects to 
the operator console 100 through a serial link 125. It is through this link 125 that 
the system control 122 receives commands from the operator which indicate the scan 
sequence that is to be performed. The pulse generator module 121 operates the system 
components to carry out the desired scan sequence. It produces data which indicates 
the timing, strength and shape of the RF pulses which are to be produced, and the 
timing of and length of the data acquisition window. The pulse generator module 121 
connects to a ?S:^....9S..3E^^;!?:}}^. amplifiers 127, to indicate the timing and shape of 
the S^aaient pulses to be produced during the ....scan . The pulse generator module 121 
also receives patient data from a physiological acquisition controller 129 that 
receives signals from a number of different sensors connected to the patient, such 
as ECG signals from electrodes or respiratory signals from a bellows. And finally, 
the pulse generator module 121 connects to a scan room interface circuit 133 which 
receives signals from various sensors associated with the condition of the patient 
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and the magnet system. It is also through the scan room interface circuit 133 that a 
patient positioning system 134 receives commands from the pulse generator module 121 
to move the patient to the sequence of desired positions to perform the scan in 
accordance with the present invention. The operator can thus control the operation 
of the patient positioning system 134 through the keyboard and control panel 102. 

Detailed Description Text (4) : 

Tne gradient waveforms produced by the pulse generator module 121 are applied to a 
gradient amplifier system 127 comprised of G.sub.x, G.sub.y and G.sub.z amplifiers. 
Each gradient amplifier excites a corresponding gradient coil in an assembly 
generally designated 13 9 to produce the magnetic field gradients used for position 
encoding acquired signals. The gradient coil assembly 139 forms part of a magnet 
assembly 141 which includes a polarizing magnet 140 and a whole -body RF coil 152. A 
transceiver module 150 in the system control 122 produces pulses which are amplified 
by an RF amplifier 151 and coupled to the RF coil 152 by a transmit/receive switch 
154. The resulting signals radiated by the excited nuclei in the patient may be 
sensed by the same RF coil 152 and coupled through the transmit/receive switch 154 
to a preamplifier 153. The amplified NMR signals are demodulated, filtered, and 
digitized in the receiver section of the transceiver 150. 

Detailed Description Text (5) : 

The transmit/receive switch 154 is controlled by a signal from the pulse generator 
module 121 to electrically connect the RF amplifier 151 to the coil 152 during the 
transmit mode and to connect the preamplifier 153 during the receive mode. The 
transmit/receive switch 154 also enables a separate RF local coil to be used during 
the receive mode. In the preferred embodiment of the invention a single, stationary 
local coil described below and shown in FIG. 7 is switched into operation. In the 
alternative, a multi -segment local coil described below and shown in FIG. 8 may also 
be used. In this case, the separate segments of the local coil are switched into 
operation by the transmit/receive switch 154 as the patient is translated to 
successive stations during the scan . 

Detailed Description Text (6) : 

The NMR signals picked up by the RF local coil are digitized by the transceiver 
module 150 and transferred to a memory module 160 in the system control 122. When 
the scan is completed and an entire array of data has been acquired in the memory 
module 160, an array processor 161 operates to Fourier transform the data into an 
array of image data. This image data is conveyed through the serial link 115 to the 
computer system 107 where it is stored in the disk memory 111. In response to 
commands received from the operator console 100, this image data may be archived on 
the tape drive 112, or it may be further processed by the image processor 106 and 
conveyed to the operator console 100 and presented on the display 104. 

Detailed Description Text (8) : 

While many pulse sequences may be used to practice the present invention, in the 
preferred embodiment a 2D gradient -recalled echo pulse sequence is used to acquire 
the NMR data. Referring particularly to FIG. 2, an RF excitation pulse 220 having a 
flip angle of 50. degree, is produced in the presence of a slab select gradient pulse 
222 to produce transverse magnetization in the 2D slab of interest (typically 100 to 
150 mm thick) . This is followed by a phase encoding gradient pulse 226 directed 
along the y axis . A readout gradient pulse 228 directed along the x axis follows and 
a partial echo (60%) NMR signal 230 is acquired and digitized as described above. 
After the acquisition, a spoiler gradient pulse 232 is applied along the z axis and 
a rewinder gradient pulse 234 is applied to rephase the magnetization before the 
pulse sequence is repeated as taught in U.S. Pat. No. 4,665,365. 

Detailed Description Text (9) : 

As is well known in the art, the pulse sequence is repeated and the phase encoding 
pulse 226 is stepped through a series of values to sample the 2D k-space in" the 
field of view, in the preferred embodiment 256 phase encodings are employed along 
the y axis. Sampling along the k.sub.x axis is performed by sampling the echo signal 
230 in trie presence of the readout 9£alierit" pulse 228 during each pulse sequence. It 
will be understood by those skilled in "the "art that only a partial sampling along 
the k.sub.x axis is performed and the missing data is computed using a homodyne 
reconstruction or by zero filling. This enables the echo time (TE) of the pulse 
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sequence to be shortened to less than 1.8 to 2.0 ms and the pulse repetition rate 
(TR) to be shortened to less than 10.0 msecs. 



Detailed Description Text (10) : 

Referring to FIG. 3, an examination of the vasculature of a patient's legs can be 
performed by dividing up the region of interest into three overlapping fields of 
view indicated at 250, 252 and 254. As shown in FIG. 4, this is accomplished by 
moving a patient table 256 to three successive locations, or stations, within the 
bore of the magnet to align the respective centers of the field of view with the 
isocenter 258 of the MR I system. A local RF coil 262 such as a 4-coil array 
consisting of anterior and posterior subunits is positioned at the system isocenter 
and it remains stationary as the patient is translated by the table 256 to three 
different stations. At each station of the table 256 the pulse sequence of FIG. 2 is 
employed to acquire NMR data from which an image of the field of view may be 
reconstructed . 



Detailed Description Text (11) : 

As shown in FIG. 7, a coil holder 200 has a closed end 202 which is attached to the 
carriage cover of the MRI system and supports an upper plate 204 and lower plate 206 
which extend into the bore of the magnet, The end of the upper plate 204 supports 
the upper RF coil unit 262A which rests on the patient and slides smoothly over the 
patient during table translation. A rectangular opening 207 in the upper plate 204 
provides space for the patient's feet. The lower plate 206 supports the lower RF 
coil unit 262B on its outer end which is positioned beneath the patient. The patient 
lies on an elevated table 208 and the lower plate 206 extends beneath the elevated 
table 208. The elevated table 208 rests on the sliding table 210 in the MRI system 
and is spaced therefrom by legs 212 to provide space for the lower plate 206. The 
patient may thus be translated between the RF coil units 262A and B which remain 
stationary at the system isocenter and in close proximity for optimal signal 
reception at all table stations. 

Detailed Description Text (12) : 

An alternative local coil structure which moves with the patient may also be used. 
Referring to FIGS. 9A and 9B, this alternative local coil is a coil array comprised 
of five coil segments 214-218, each comprised of four coil elements. The coil 
segments 214-218 are positioned on the patient and distributed along the entire 
region of interest to be imaged - - in this case, the legs and feet. The coil elements 
are supported by fabric (not shown) which is sewn into a pant-like garment that 
clothes the patient's legs. 

Detailed Description Text (13) : 

Each coil segment 214-218 has four coil elements that acquire NMR data from the FOV . 
Two of the coil elements are positioned on top, or anterior, of the lower extremity 
as seen in FIG. 9B, and the other two elements are positioned below, or posterior, 
of the lower extremity. The coil elements connect together to form one of the coil 
segments 214-218, and each is separately connected through terminals 219 to the 
above-described transmit/receive switch 154. Coupling between coil segments 214-218 
is of little concern, because only one segment is operative at any moment during the 
scan . As the patient is translated to bring the successive coil segments 214-218 
into the FOV of the MRI system, the aligned coil segment is connected to the system 
transceiver module 150. MRA data is acquired using that coil segment and the patient 
is then translated to the next station and the next coil segment is enabled. 

Detailed Description Text (14) : 

The multi -element coil segments may be constructed using a number of well-known 
methods. Coupling between the separate coil elements is minimized to increase the 
SNR of the data acquired using the coil segment. Such decoupling is achieved using 
the methods taught in U.S. Pat. No. 4,825,162 issued to Roemer, et al . , and entitled 
"Nuclear Magnetic Resonance (NMR) ] : w®3jj}3.. with Multiple Surface Coils" or U.S. Pat. 
No. 4,721,913 issued to Hyde, et al .""and "entitled "NMR Local Coil Network". 

Detailed Description Text (15) : 

As shown in FIG. 6, the scan i s performed under the dir^ection of a stored program 
which directs the pulse generator module 121 to carry out a sequence of stegs • As 
indicated by pSpcess block 370, the pulse control module 121 sends command signals 
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to the patient positioning system 134 which moves the table to align the first FOV 
at the system isocenter 258. A loop is then entered in which NMR data is acquired 
for a complete image . More specifically, the pulse control module 121 directs the 
MR I system to perform the pulse sequence of FIG. 2 to acquire one view of NMR data 
as indicated at process block 372. Phase encodings are stepped as described above 
and additional views are acquired until all of k- space has been sampled and the 
acquisition is complete as determined at decision block 374. 

Detailed Description Text (16) : 

At the completion of the first image acquisition the pulse control module 121 
commands the patient positioning system 134 to move the table 256 to the next 
station to align the next field of view at the system isocenter 258 as indicated at 
process block 376. In the preferred embodiment this distance is 25 cm. Another image 
is then acquired at 372 and 374. This sequence of acquiring an image and moving the 
patient table 256 continues until the last field of view in the scan has been 
acquired as determined at decision block 378. Both the FOV sizes and the separation 
between the anterior and posterior coils may vary from FOV to FOV for optimal image 
resolution and signal-to-noise ratio. 

Detailed Description Text (17) : 

The acquired field of view images are registered with each other and combined to 
form a single image of the much larger region of interest. The patient is 
immobilized to the table 256 with straps (not shown) to minimize misregistration 
between the field of view images . In addition, a marker made of a contrast agent 
solution which produces a high level signal can be placed along side the patient and 
used to align each field of view image such that they are in registration with each 
other and form a single, contiguous image of the entire region of interest. 

Detailed Description Text (18) : 

In the preferred method the acquisition of each FOV is timed to occur when maximum 
image contrast is produced by a contrast agent administered intravenously to the 
patient. A test bolus (eg. 5 ml of gadolinium) is administered intravenously to 
measure the bolus timing curve at a predetermined region of interest, typically the 
second FOV . As shown in FIG. 5, for example, the curve 274 peaks in the first FOV 
250 at approximately 23 seconds after the contrast agent is administered, the curve 
275 peaks in the second FOV 252 at approximately 29 seconds after the contrast agent 
is administered and the curve 276 peaks in the third FOV 254 at approximately 35 
seconds after contrast administration. This establishes when the data should be 
acquired for best image contrast. 

Detailed Description Text (19) : 

If the three FOVs 250, 252, 254 are acquired in order, the intraarterial column of 
contrast agent is imaged as it transits through the patient's legs. To maximize 
image contrast, the velocity of the bolus movement is estimated and the table 
translation is timed to acquire each FOV at peak contrast. One way to estimate the 
bolus velocity is to measure the arterial velocity (averaged over the cardiac cycle) 
by acquiring a phase contrast image through the artery of interest in the center of 
the region of interest as described in U.S. Pat. No. Re 32,701. Blood flow velocity 
can be determined from this image and an estimate made of the bolus velocity. 
Another way to estimate the bolus velocity is to perform a timing scan in which the 
bolus advancement through the FOV is monitored using a high temporal frame rate. 
This second method provides a rapid estimation of bolus velocity that is accurate 
enough to use with the present invention. From this velocity estimate and the peak 
contrast timing for the central FOV 252, the timing for the acquisition of the FOVs 
250 and 254 can be calculated. These timings will determine the rate of table 
translation between FOV acquisitions. 

Detailed Description Text (20) : 

The timing for table translation is estimated from three parameters: the estimate of 
bolus velocity; the duration of data acquisition at each j?0V; and the time for 
moving the patient from one FOV to the next. Optimal results are achieved in the 
preferred embodiment if the leading edge of the bolus contrast agent has flowed into 
the next FOV for about 5 sec (T . sub . 12*) ." The time for each additional FOV (T.sub.s) 
is the time for table translation from one FOV to the next (3 seconds or less if 
performed manually) , and the Af.^SArtS. duration (about 5-7 sec) . The ..total, length of 
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the contrast agent bolus column (L.sub.b) required to perform the scan can then be 
estimated from the following equation: 

Detailed Description Text (21) : 

Here N.sub.s is the number of FOVs to be imaged, and V.sub.b is the estimated bolus 
velocity. L.sub.s is the table translation distance between stations, which is 
slightly smaller than FOV to provide overlapping acquisitions between adjacent 
stations. A constant bolus velocity is used in the above equation. This concept can 
also be extended straight forwardly to the case of varying bolus velocity using 
integration. The imaging duration at each station is kept at several seconds to 
account for velocity variation as the contrast column travels to the more peripheral 
FOVs . Accordingly, the contrast duration is T.sub.b =L.sub.b /V.sub.b and the 
injection duration is less than this by the cardiac dispersion effects (.tau..sub.c 
.about. 5 seconds). For a given injection rate, (R.sub.i) the required contrast dose 
is as follows: 

Detailed Description Text (22) : 

An example of using this timing and contrast dose information is illustrated in FIG. 
8. In this example, the length of each imaging region is 25 cm and indicated by 
shaded region 290, and it is translated to the second, third and fourth FOVs at 30, 
40 and 50 seconds respectively. The contrast agent bolus 292 transits the region at 
a speed of 3 cm/sec and imaging of the first FOV continues until the leading edge of 
the contrast agent bolus 2 92 is 5 seconds into the second FOV . The time at each 
additional FOV is 10 seconds, 7 seconds for acquiring image data and 3 seconds for 
moving the next FOV into the imaging region 290. The total length of the contrast 
agent bolus is 55 cm, which corresponds to a bolus duration of 18.3 seconds. A bolus 
injection duration of 15 seconds is used, and at an injection rate of 1 cc/sec, the 
total dose is 15 cc (using 3.3 sec . tau . . sub . c) . 



Detailed Description Text (23) : 

As is well known in the art, "mask" images may also be acquired before contrast 
agent is administered. In the preferred embodiment, the mask images are acquired 
before bolus injection at each of the stations. The above described contrast 
enhanced acquisition is then performed and a complex subtraction of corresponding 
voxels in the mask image from the reconstructed image is performed to further 
enhance image quality. 

Detailed Description Text (24) : 

The preferred method for practicing the invention includes conducting a "scout scan " 
to confirm that the vasculature of interest is within the FOV of at least one of the 
stations. This is followed by a "timing scan " in which information is acquired 
regarding the arrival and velocity of the contrast agent bolus, and by a "mask scan " 
in which the mask image data is acquired. And finally, the "bolus tracking MRA scan " 
is conducted as described above and shown in FIG. 6. The following is an exemplary 
procedure for carrying out these steps in a four station scan of a patient's legs. 

Detailed Description Text (25) : 
Scout Scan 

Detailed Description Text (26) : 

A multiple location, sagittal acquisition using a 2D fast gradient echo sequence 
( FOV =40 cm, TE/TR=2/9, receiver bandwidth=16 kHz, 256 . times . 256 , 1/2 FOV , total 
acquisition time .about. 15 sec) is performed at the first station. The patient is 
translated to the second station and the same scan is repeated immediately without 
any prescan. This process repeats until all stations are imaged . For a 4 -station 
scout scctn, the total imaging time is about 70 sec (allowing 10 seconds for each 
table motion) . 

Detailed Description Text (27) : 
Timing Scan 

Detailed Description Text (28) : 

A thick coronal slab timing acquisition is typically performed at the second station 
using a gadolinium contrast agent test dose (3-5 mL, at 1 mL/sec injection rate) and 
a 2D fast .gradient^ echo sequence { FOV=32-40 cm, 130-170 mm thick, TE/TR=2/9, 
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receiver bandwidth=16 kHz, 256 . times . 160-192 , 40-60 image frames, total acquisition 
time .about. 60-90 sec to ensure the capture of the contrast peak). From this timing 
acquisition, the arrival time of the contrast agent in the FOV is measured, and the 
bolus velocity is estimated based on the time required for the bolus to transit the 
FOV . 

Detailed Description Text (29) : 
Mask Scar: 

Detailed Description Text (30) : 

Similar to the timing scan, a thick coronal slab containing all arteries of interest 
is acquired using a fast 2D gradient echo sequence ( FOV =32 cm, 130-170 mm thick, 
TE/TR=2/9, receiver bandwidth=16 kHz, 256 . times . 160-192, 40-60 image frames, total 
acquisition time .about. 60-90 sec. In contrast to the timing scan, the table is 
translated to each station at times determined by the timing scan and mask data is 
acquired at each FOV . 

Detailed Description Text (31) : 
Bolus Tracking MRA Scan 

Detailed Description Text (32) : 

The contrast agent is administered and a scan is conducted using the same 
prescription as the mask scan . As with the mask scan, the data acquisition is begun 
at the time, and the table is translated at the times determined by the timing scan . 



Detailed Description Text (34) : 

Next, a multiple station acquisition without contrast injection was performed to 
obtain mask data. Immediately following the mask acquisition, a gadolinium-enhanced 
bolus tracking acquisition was performed with a gadolinium dose of 0.1-0.2 mmol/kg. 
The table translation from the first station to the second station occurred at the 
peak time of the contrast curve at the distal location of the second station. Table 
translations were 20-30 cm per step, providing 2-12 cm overlap. Hand injection and 
later an MR compatible power injector (Spectris MR Injector, MedRad, Pittsburgh, 
Pa.) was used. All infusions were followed immediately by a 30 mL saline flush. 

Detailed Description Text (35) : 

In one patient, for example, the timing bolus arrived at the popliteal artery 
approximately 27 sec after injection, and it took approximately 6 sec to traverse 
the central station. Three-station acquisitions were planned for this subject. A 15 
mL dose at 1 mL/sec injection rate was used for bolus tracking acquisition (15 sec 
injection time) . Accordingly, table translation from the first station to the second 
station occurred between 34-36 sec after injection, allowing about 10 sec for 
imaging the first station. Imaging time at the second station was 12 sec. Then the 
table was translated to the third station between 48-50 sec after injection. The 
total acquisition time was 60 sec. All arteries from the external iliac artery in 
the pelvis through the anterior tibial and posterior tibial arteries at the ankle 
are well depicted. 

Detailed Description Text (36) : 

Our results support the hypothesis that MRA of the lower extremity can be performed 
rapidly using a bolus tracking strategy, a stepping table, and a 2D MR imaging 
technique. With this bolus tracking MR technique, MR angiography of the entire lower 
extremity can be performed in approximately one minute. Including time for acquiring 
the scout, bolus timing and mask data sets, the total examination time for MRA of 
the entire lower extremity is between ten and fifteen minutes. This is a substantial 
reduction in scan time from that of the standard TOF technique (typically 60-90 
minutes) . 

Detailed Description Text (37) : 

The 2D acquisition employed in this study provides high temporal resolution while 
complex subtraction of the mask injac^e provides high contrast-to-noise ratio with 
minimal artifacts . A frame rate of 1-2 sec per image can be maintained while 
providing 1 mm. times. 1 mm in-plane resolution with good SNR. This rapid acquisition 
enables iiTia^ging^ of the same contrast bolus sequentially as it transits the arterial 
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tree * Xtnaginq of each station for five to ten seconds permits the acquisition of the 
contrast enhanced lower extremity arteries in 30 seconds. This is less than the time 
required for imaging one station using a 3D acquisition. Consequently, a much 
smaller contrast dose can be used in bolus tracking 2D MR. An alternative way to 
take advantage of 2D acquisition is to acquire multiple projections at each station 
by alternating slab orientation. Because the in-plane resolution is higher than that 
of a MIP image from a 3D acquisition, this may be a more effective way to obtain 
additional spatial information of the targeted arteries. 

Detailed Description Text (38) : 

An alternative method to the "bolus tracking" of the preferred embodiment is to 
separately administer contrast agent prior to each FOV acquisition. In this case a 
mask image is acquired just prior to each bolus injection, the bolus is 
administered, and the image acquisition for the particular FOV is performed when 
peak contrast is obtained in that FOV . The table is then translated to the next FOV 
ana the process is repeated. 

Detailed Description Text (39) : 

Other data processing techniques may be used to generate angiograms . Methods other 
than complex subtraction, such as digital filters, may be used to generate 
angiograms. The criteria for digital filter design is to extract dynamic contrast 
(or signal variation) from a temporal series of images f and vessel contrast is 
determined by this dynamic contrast. For example, matched filters can be used to 
generate an angiogram that is a sum of individual arteriograms generated by complex 
subtraction . 

Detailed Description Text (40) : 

Another application for the invention is to image the aorta in the body trunk. The 
aortic arch and the abdominal aorta can be image in two separate FOVs using a fast 
3D acquisition technique (temporal resolution<one 3D volume per 20 sec) and contrast 
infusion. In addition, the table stepping between different FOVs can be an important 
technical component for general fluoroscopic MR imaging and MRI -guided 
interventions . 

CLAIMS : 

1. A method for producing an image with an MRI system, the steps comprising: 

a) positioning a patient in the MRI system; 

b) acquiring NMR image data from the patient over a first field of view while the 
patient is stationary relative to the MRI system; 

c) translating the patient in the MRI system; 

d) acquiring NMR image data from the patient over a second field of view while the 
patient is stationary relative to the MRI system; 

e) registering the NMR image data acquired over the first and second field of views 
and reconstructing an image over a region of interest which includes the first and 
second field of views. 

2 . The method as recited in claim 1 in which a contrast agent is administered to the 
patient and the acquisition of NMR image data over the first and second field of 
views occurs when the contrast agent provides substantially maximum contrast 
therein . 

5. The method as recited in claim 1 in which the patient is translated additional 
times, NMR image data is acquired from additional field of views, and the 
reconstructed image includes NMR image data from the additional field of views. 

6. The method as recited in claim 1 in which the patient is translated a preselected 
distance and the NMR image data acquired over the first field of view is registered 
with the NMR image data acquired over the second field of view by establishing their 
relative positions using said preselected distance. 
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7. The method as recited in claim 1 in which the NMR image data is acquired using a 
local RF coil which is positioned adjacent the patient and which remains stationary 
as the patient is translated. 

8. The method as recited in claim 1 in which the NMR image data is acquired using a 
local RF coil having a plurality of coil segments which are supported by the 
patient, and as the patient is translated different coil segments are switched into 
operation to acquire NMR image data. 

9. The method as recited in claim 1 in which a contrast agent is administered to the 
patient and a timing scan is conducted to estimate the optimal time at which step b) 
is to be performed, and to estimate the velocity of the contrast agent as it flows 
through one of said field of view. 

10. The method as recited in claim 9 in which the performance of step c) is 
determined in part by the estimated velocity of the contrast agent. 

11. A method for producing an image with an MR I system, the steps comprising: 

a) performing a timing scan in which a contrast agent is administered to a patient 
and the arrival time and the velocity of the contrast agent in a region of interest 
is estimated from NMR data acquired by the MR I system; 

b) performing a bolus tracking scan to acquire NMR data from the region of interest 
by translating the patient to a plurality of stations and acquiring NMR image data 
from a corresponding plurality of field of views within the region of interest; and 

c) reconstructing the image from NMR data acquired from all of said field of views; 

wherein the start of the bolus tracking scan is determined by said estimated arrival 
time and the translation of the patient to said plurality of stations is timed in 
part by said estimated velocity. 
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